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Abstract 

Durability  and  dynamic  failure  properties  are  critical  parameters  for  naval  composite 
ships  in  seawater.  However,  previous  measurements  have  significantly  underestimated  the 
intrinsic  durability  of  composite  structures.  In  this  study,  a  new  composite  “fish  tank"  approach 
was  developed  to  accurately  simulate  the  material/mechanics  conditions  of  composite  structures 
in  seawater.  Four  E-glass/vinyl  ester  composite  specimens  were  weakly  bonded  together  to  form 
a  seawater  tank.  Therefore,  only  one  surface  of  each  composite  specimen  was  exposed  to 
seawater  and  this  surface  was  subjected  to  drop-weight  impact,  which  is  very  similar  to  the 
dynamic  failure  process  of  ship  structures  subjected  to  underwater  explosion.  A  series  of 
experiments  of  dry  and  wet  composite  specimens  were  conducted  to  characterize  the  impact 
damage,  and  the  residual  compression  strength.  Results  of  more  than  50  specimens  with  more 
than  21 -month  seawater  exposure  show  that  the  compression  after  impact  strengths  of  the  wet 
specimens  reduced  less  than  9%  compared  to  the  dry  baseline  specimens.  Therefore,  durability 
of  current  marine  composite  is  much  better  than  previously  measured.  Our  simple  and  novel 
“fish  tank”  durability  experiment  is  expected  to  provide  accurate  justification  for  the  required 
maintenance  period  of  new  composite  ship  structures. 

Technical  Section 

Technical  Objectives 

The  main  objectives  of  this  project  include:  1)  developing  a  novel  approach  to  accurately 
simulate  the  material/mechanics  conditions  of  composite  structures  in  seawater;  2)  characterizing 
the  impact  damage,  and  the  residual  compression  strength  of  marine  composites  as  a  function  of 
time  of  seawater  exposure;  and  3)  conducting  a  combined  experimental  and  numerical 
investigation  of  the  interface  shear  strength  measurements. 


Technical  Approach 

In  this  study,  a  new  composite  “fish  tank”  approach  was  developed.  Four  E-glass/vinyl 
ester  composite  specimens  were  weakly  bonded  together  to  form  a  seawater  tank.  Therefore, 
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only  one  surface  of  each  composite  specimen  was  exposed  to  seawater  and  this  surface  was 
subjected  to  drop-weight  impact,  which  is  very  similar  to  the  dynamic  failure  process  of  ship 
structures  subjected  to  underwater  explosion. 

Materials  and  sample  preparation 

Glass  fiber  reinforced  vinyl  ester  (glass/VE)  panels  were  produced  using  vacuum  assisted 
resin  transfer  molding  (VARTM)  by  Prof.  U.  Vaidya's  group  at  the  University  of  Alabama  at 
Birmingham.  Eight  layers  of  plain  weave  glass  fabric  (CWR  2400/50  plain  weave.  Composites 
One,  LLC)  were  used  to  produce  the  panels  with  approximately  5mm  thickness  which  is  required 
by  ASTM  D  7137  “Standard  Test  Method  for  Compressive  Residual  Strength  Properties  of 
Damaged  Polymer  Matrix  Composite  Plates".  The  fiber  fraction  of  the  panels  was  found  to  be 
54%  vol.  after  bum  off  testing  was  conducted.  Compression  after  impact  (CAI)  testing  samples 
with  a  dimension  of  101.6  mm  x  152.4  mm  (4"  x  6")  were  cut  and  machined  to  meet  the  strict 
dimension  requirement  specified  in  ASTM  D  7137. 


Construction  of  a  composite  seawater  tank 


Figure  1 .  A  composite  tank  before  construction  (left)  and  after  construction  with  seawater  inside  (right) 


As  shown  in  Figure  1,  silicone  rubber  as  aquarium  sealant  (Perfecto  Manufacturing, 
Noblesville,  IN)  were  applied  to  four  slots  of  a  base  PMMA  plate  before  four  composite 
specimens  were  inserted.  PMMA  has  very  little  reaction  with  seawater.  The  reason  to  use 
silicone  rubber  is  that  it  provides  enough  bonding  strength  under  water  pressure.  Meanwhile,  it  is 
not  too  strong  for  to  break/debond  this  tank  into  four  composite  plates  for  future  impact 
experiments.  After  one  week  of  the  construction  of  this  tank  (full  bonding  strength),  it  was  filled 
with  synthetic  seawater  (Ricca  Chemical  Co.,  TX).  This  tank  will  be  disassembled  after  certain 
periods  of  time  such  as  three  months,  six  months  etc.  to  conduct  impact  and  compression 
experiments  (see  Figure  2).  The  impact  experiments  of  dry  specimens  were  conducted  to  provide 
baseline  data. 
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Figure  2.  Layered  composite  specimens  subjected  to  out-of-plane  impact  and  compression 


Figure  3.  Typical  impact  damage  on  the  front  and  back  surfaces,  (a)  and  (b), 
and  typical  compression  failure  of  the  impacted  specimen  (c)  and  (d). 

Impact  experiments  of  durability  specimens 

Impact  damage  was  introduced  using  a  drop  tower  setup.  All  samples  (fixed  four  edges) 
were  subjected  to  an  impact  (60  joules  impact  energy)  using  a  16mm  (5/8")  diameter  hemisphere 
impactor.  Damage  zones  of  the  impacted  samples  are  clearly  seen  in  Figure  3(a),  (b).  For  the 
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front  surface  directly  subjected  to  impact,  dark  areas  represent  internal  delamination,  with 
possible  several  delaminations  at  the  different  interfaces.  As  discussed  by  Xu  and  Rosakis  [1], 
these  delaminations  are  mainly  shear-dominated,  so  the  interlaminar  shear  strength  is  an 
important  parameter  for  delamination  resistance  characterization.  Also,  two  major  matrix  cracks 
were  observed  near  the  impact  site  (as  shown  by  two  dark  mark  lines).  One  matrix  crack  was 
along  the  horizontal  direction  and  the  other  one  was  along  the  vertical  direction.  On  the  back 
surface  of  the  impacted  specimen,  fiber  breakage  was  observed  at  the  impact  site  and  this  failure 
mode  contributed  to  major  impact  energy  absorption.  Meanwhile,  fiber/matrix  debonding 
appeared  as  white  thin  lines  on  the  back  surface  of  the  impacted  specimens.  These  four  major 
failure  modes  indeed  make  different  contributions  to  the  composite  impact  resistance  [2],  and  we 
believe  fiber  breakage  and  delamination  play  the  major  role  to  absorb  impact  energy. 

Compression  tests  for  impacted  specimens 


70000 
60000 

~  50000 

z 

•p  40000 

es 

-2  30000 

20000 
10000 
0 

0.0  0.5  1.0  1.5  2.0 

Displacement  (mm) 

Figure  4.  A  ty  pical  load-displacement  curve  of  a  compressive  test  of  an 
impacted  marine  composite  laminate 

Impacted  samples  were  mounted  into  a  compression  fixture.  A  loading  rate  of  1  mm/min  was 
used.  The  progressive  compression  failure  started  from  the  impact  damage  as  shown  in  Figure  4. 
Initially,  as  the  compression  load  increased,  delamination  from  the  previous  impact  propagated 
in  a  local  buckling  form.  Unlike  impact-induced  delamination,  its  propagation  is  mainly  opening- 
dominated.  Delamination  also  appeared  along  the  horizontal  matrix  crack  and  this  matrix  crack 
extended  to  the  two  edges  as  the  compressive  loading  increased,  as  seen  in  Figure  3.  The  final 
failure  (maximum  load)  was  controlled  by  a  shear  crack  near  the  horizontal  matrix  crack  as  seen 
in  Figure  4.  An  inclined  angle  around  30-45  degrees  (with  respect  to  the  compressive  loading 
direction)  was  observed  from  the  two  vertical  edges  of  the  failed  specimen.  These  results  are 
similar  to  previous  compressive  failure  results  reported  by  Daniel  [3],  Tsai  and  Sun  [4],  Oguni 
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and  Ravichandran  [5].  A  load-displacement  curve  is  illustrated  in  Figure  4  for  a  compressive 
experiment  of  an  impacted  specimen.  The  initial  non-linear  part  is  caused  by  the  initial  gap  of  the 
compressive  fixture.  Then  a  long  linear  load-displacement  part  is  recorded.  The  failure  mode 
starts  from  the  opening  delamination  from  the  impacted-induced  delamination  (shear- 
dominated),  followed  by  a  sudden  propagation  of  the  longitudinal  matrix  crack  and  a  final  shear 
crack  appeared  along  the  specimen  edge  based  on  the  recorded  high-definition  video. 

A  master  curve  of  the  accurate  durability  measurement  using  our  new  approach 

Table  1  shows  all  CAI  data  up  to  21  months.  The  dry  specimen  is  used  as  a  baseline 
specimen  or  comparison.  Since  the  CAI  strength  combines  the  effects  of  the  seawater  exposure 
and  impact  damage,  it  is  very  convenient  to  use  this  simple  parameter  to  present  complicated 
impact  damage  and  the  seawater  effect. 


Table  1.  History  of  the  Compression- After-Impact  (CAI)  Strength  in  seawater 


Time  (months) 

Mean  CAI  (Mpa) 

%  Reduction 

0 

132.98  ±7.59 

0 

4 

140.40  ±4.37 

+5.58 

9.5 

121.39  ±  10.98 

-8.72 

13 

125.63  ±9.05 

-5.53 

17 

128.15  ±8.77 

-3.63 

21 

122.22  ±9.22 

-8.09 

From  the  table,  we  notice  that  the  CAI  strength  reduction  is  less  than  9%  after  one-year 
seawater  exposure.  This  is  much  less  than  40%  CAI  strength  reduction  as  reported  by  Imielinska 
and  Guillaumat  [6]  on  the  same  compression  after  impact  experiments  with  different  composite 
materials.  Obviously,  if  the  durability  of  a  new  composite  ship  reduced  so  dramatically  in  less 
than  one  year,  this  ship  should  not  be  manufactured.  This  comparison  confirms  that  our  new 
approach  produces  more  reasonable  data  as  our  experiments  simulate  the  right  material 
conditions. 

These  CAI  data  are  also  plotted  in  Figure  5.  A  slight  increase  in  CAI  strength  for  the 
specimen  after  four-month  seawater  exposure  is  probably  due  to  the  specimen  size  effect.  The 
average  thickness  of  this  set  of  specimens  is  at  least  10%  higher  than  other  specimens.  The  CAI 
strength  is  not  a  material  property'  so  it  is  sensitive  to  the  specimen  sizes  especially  the  specimen 
thickness.  After  four-month  seawater  exposure,  the  CAI  strength  reduction  has  very  small 
change  up  to  21  months,  which  might  be  related  to  seawater  saturation  inside  the  specimens. 
Since  no  additional  seawater  enters  the  composite  specimens  after  a  certain  amount  of  time  [7], 
no  new  damage  caused  by  seawater  effect  such  as  fiber-matrix  interface  strength  reduction  will 
occur.  Therefore,  the  CAI  strength  of  the  wet  specimen,  which  combines  the  effects  of  impact 
damage  (not  a  function  of  seawater  exposure  time),  and  the  seawater  damage  (such  as  to  the 
matrix),  will  have  very  small  change  in  one  to  two  years.  Since  there  is  no  protective  coating  or 
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painting  on  the  composite  plates,  the  future  composite  ship  structures,  which  have  strong  coating 
or  painting,  should  have  excellent  durability  and  dynamic  failure  properties. 
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Figure  5.  A  master  curve  of  the  compressive  strength  change  as  a  function  of 
seawater  exposure  time 


Experimental  studies  of  two  common  approaches  of  shear  strength  measurements 

Since  the  previous  impact  failure  modes  included  shear-dominated  delamination  and 
shear  cracks,  shear  strength  measurements  are  very  important  for  composite  materials.  Among 
several  shear  strength  approaches,  the  losipescu  shear  test  has  been  used  for  the  past  20-30  years 
because  the  losipescu  shear  specimen  provides  a  test  region  of  uniform  shear  stress,  rather  than  a 
parabolic  variation  shown  by  the  short-beam  shear  specimen  as  seen  in  Figure  6. 
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Figure  6.  Schematic  of  (a)  Butt-joint  shear  specimen  and  (b)  Iosipescu  shear  specimen. 


However,  due  to  different  shear  stress  distributions  of  these  two  approaches,  will  they  lead  to 
different  shear  strength  measurements?  No  one  answers  this  fundamental  question  and  we  try  to 
provide  partial  answer  using  simplified  materials,  i.  e.,  bonded  homogenous  materials  with  weak 
interfaces.  Test  specimens  were  made  of  three  different  types  of  materials  including  Aluminum, 
PMMA  and  Polycarbonate.  The  reason  for  using  bonded  specimens  is  that  we  want  the  shear 
specimens  to  fail  at  the  interface  by  shear  and  not  in  tension  at  the  upper  edges.  The  adhesive 
used  in  this  experiment  is  Loctite  384.  One  of  the  key  requirements  in  choosing  the  adhesive  was 
that  its  strength  should  be  close  to  or  lesser  than  that  of  the  bulk  material.  The  individual  halves 
were  machined  to  ensure  that  they  were  of  the  desired  dimensions  and  were  sand-blasted  at  the 
bonding  surfaces  to  provide  good  adhesion.  The  specimens  were  bonded  using  a  fixture  to 
guarantee  their  dimensionality  after  bonding  under  room  temperature. 

The  test  set-up  consists  of  three  parts  including  a  mechanical  system  for  testing  the 
specimens,  an  optical  system  used  to  develop  photoelasticity  fringe  patterns  and  an  imaging 
system  to  record  and  capture  the  images.  All  specimens  were  loaded  to  failure  and  their  strength 
data  were  recorded  as  shown  in  Table  2.  It  should  be  noted  that  the  failure  of  all  specimens  was 
sudden  and  brittle.  From  the  above  table  it  is  seen  that  the  average  values  are  close  and  the  data 
from  the  Iosipescu  shear  test  are  slight  higher.  The  Iosipescu  specimen  shows  a  higher  value  of 
measured  bonding  strength  because  of  the  lower  cross  sectional  area  involved.  The  Iosipescu 
specimen  has  about  60%  of  the  cross  sectional  area  of  the  butt-joint  shear  specimen  and  thus  it 
leads  to  a  greater  value  of  the  measured  bonding  strength  due  the  “size  effect”.  These  results 
show  that  the  shear  stress  distribution  has  least  effect  on  the  shear  strength  measurements.  Of 
course,  this  result  is  only  applicable  to  any  interface  without  stress  singularity  such  as  caused  by 
a  crack  (macro-scale)  or  bi-material  comer. 
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Table  2.  Measured  Bonding  Shear  Strength  Data  fat  least  20  specimens  tested  for  each  case) 


Materials 

Iosipescu  shear  (MPa) 

Butt-joint  shear  (MPa) 

Difference  (%) 

Aluminum-Al 

10.75 

10.16 

+  5.8 

Polycarbonate-PC 

10.99 

8.51 

+  23.3 

PMMA-PMMA 

11.58 

10.19 

+  13.6 

Modeling  and  simulation  of  the  shear  strength  comparison 

A  finite  element  model  of  the  specimen  was  built  using  ANSYS  11.0  to  obtain  the 
distribution  of  the  stresses  in  the  specimen.  Two  different  types  of  materials  were  chosen  to 
represent  the  actual  experiment.  A  linear,  elastic  and  isotropic  material  was  chosen  as  only  brittle 
materials  are  considered  in  this  analysis.  The  specimen  was  meshed  using  Plane  42  quad 
elements.  A  finer  mesh  was  used  in  the  notch  region  in  order  to  capture  any  stress  singularities. 
A  coarser  mesh  was  used  elsewhere  to  save  on  computational  time  and  effort.  Plane  42  elements 
were  exclusively  used  in  order  to  meet  to  meet  the  requirements  of  the  plotting  software  Tecplot. 
In  order  to  incorporate  a  realistic  simulation  of  the  load  transfer  from  the  fixture  to  the  specimen, 
an  iterative  procedure  was  adopted.  The  left  portion  of  the  specimen  is  fixed  and  cannot  move  in 
the  vertical  dimension.  Displacement  constraints  are  applied  similar  to  what  is  observed  in 
reality.  The  top  left  edge  node  is  restrained  from  moving  in  both  x  and  y  directions  in  order  to 
prevent  rigid  body  movements.  On  the  right  portion,  an  initial  vertical  uniform  displacement  is 
applied.  Now  the  reaction  loads  at  each  of  the  nodes  with  a  constraint  is  checked  to  verify  that 
they  are  not  in  tension.  The  constraint  was  removed  from  those  nodes  which  showed  a  tensile 
reaction  force  instead  of  a  compressive  one  and  a  new  analysis  with  the  updated  set  of  boundary 
conditions  was  carried  out.  This  was  done  until  convergence  and  the  total  applied  load  was 
obtained  by  integrating  the  shear  stress  across  the  interface  of  the  specimen. 

The  load  applied  in  the  form  of  displacements  was  initially  obtained  from  the 
experimental  load-displacement  plot.  The  actual  load  on  the  specimen  was  later  obtained  by 
integrating  the  shear  stress  along  the  interface.  This  procedure  was  iteratively  carried  until  the 
required  load  (for  matching  the  fringe  patterns)  was  obtained.  The  numerical  modeling  of  the 
shear  specimens  w  as  validated  by  comparison  of  the  numerical  photoelastic  fringe  plots  with  the 
experimental  ones.  The  comparison  for  an  Iosipescu  shear  specimen  is  shown  in  Figure  7.  It  is 
seen  that  the  experimental  and  numerical  patterns  have  a  pretty  close  match.  The  fringes  appear 
concentrated  around  the  top  and  bottom  loading  blocks.  This  is  due  to  the  highly  non-linear  force 
distribution  on  the  upper  right  and  lower  left  loading  blocks. 
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These  results  show  that  the  shear  stress  distribution  has  least  influence  on  the  strength 
measurements.  Therefore,  future  simplified  and  efficient  methods  can  be  developed  to  measure 
important  shear  strengths  as  long  as  the  shear  stress  is  not  singular. 


(c)  Applied  Load  =  200  +/-  5  N  (d) 

Figure  7.  Comparison  between  experimental  and  finite  element  generated  fringe 
patterns  for  an  Iosipescu  specimen  made  of  bonded  Polycarbonate 
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Real-Time  Experimental  Investigation 
on  Dynamic  Failure  of  Sandwich  Structures 
and  Layered  Materials 


L.  Roy  Xu  and  Ares  J.  Rosakis 


Abstract  We  present  a  systematic  experimental  investigation  of  the  generation  and 
subsequent  evolution  of  dynamic  failure  modes  in  sandwich  structures  and  layered 
materials  subjected  to  out-of-plane  low-speed  impact.  Model  sandwich  specimens 
involving  a  compliant  polymer  core  sandwiched  between  two  metal  layers  and  other 
model  layered  materials  were  designed  to  simulate  failure  evolution  mechanisms  in 
real  sandwich  structures  and  layered  materials.  High-speed  photography  and  dy¬ 
namic  photoelasticity  were  utilized  to  study  the  nature  and  sequence  of  such  failure 
modes.  In  all  cases,  inter-layer  (interfacial)  cracks  appeared  first.  These  cracks  were 
shear-dominated  and  w-ere  often  intersonic  even  under  moderate  impact  speeds. 
The  transition  from  inter-layer  crack  growth  to  intra-layer  crack  formation  was 
also  observed.  The  shear  inter-layer  cracks  kinked  into  the  core  la>er,  propagated 
as  opening-dominated  intra-layer  cracks  and  eventually  branched  as  they  attained 
high  enough  growth  speeds  causing  brittle  core  fragmentation.  In-depth  failure  me¬ 
chanics  experiments  on  the  dynamic  crack  branching,  crack  kinking  and  penetration 
at  a  weak  interface,  interfacial  debonding  ahead  of  a  main  incident  crack  were  also 
conducted  to  understand  the  physical  insight  of  the  dynamic  failure  modes  and  their 
transition  observed  from  sandwich  structures  and  layered  materials. 


1  Introduction 

Layered  materials  and  sandwich  structures  have  diverse  and  technologically  inter¬ 
esting  applications  in  many  areas  of  engineering.  These  include  the  increased  use 
of  composite  laminates  in  aerospace  and  automotive  engineering;  the  introduction 
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of  layered  concrete  pavements  in  civil  engineering;  the  use  of  thin  films  and  lay¬ 
ered  structures  in  micro-electronic  components,  and  recently,  the  introduction  of 
sandwich  structures  in  a  variety  of  naval  engineering  applications  (Hutchinson  and 
Suo  [1];  Rajapakse  [2]).  While  failure  characteristics  of  layered  materials  and  sand¬ 
wich  structures  subjected  to  static  loading  have  been  investigated  extensively  in  the 
past  years,  their  dynamic  counterparts  have  remained  elusive  (Sun  and  Rechak  [3]; 
Cantwell  and  Morton  [4];  Abrate  [5]).  Indeed,  the  presence  of  highly  complex  and 
transient  dynamic  failure  modes  in  such  materials  and  the  inaccessibility  of  inter¬ 
nal  damage  have  resulted  in  experimental  studies  limited  to  only  the  final  impact 
damage  and  the  residual  strengths.  To  begin  addressing  the  need  for  real-time  ob¬ 
servations  of  failure  events,  the  work  presented  here  focuses  on  the  study  of  such 
events  in  model  sandwich  structures,  and  in  particular,  on  the  identification  of  their 
nature,  chronological  evolution  and  interaction. 

To  identify  the  evolution  of  failure  modes,  it  is  convenient  to  first  classify 
these  modes  based  on  the  material  constitutions  of  layered/reinforced  structures. 
As  shown  schematically  in  Fig.  1 ,  there  are  two  major  categories  of  failure  observed 
in  post-mortem  studies.  The  first  major  failure  category  is  decohesion  (or  cracking) 
between  bonded  layers  at  an  interface.  This  is  often  referred  to  as  delamination  in 
composite  laminates  or  interfacial  debonding  in  thin  films  or  sandwich  structures.  It 
is  also  called  inter-layer  failure.  Generally,  two  distinct  inter-layer  failure  modes  are 
observed.  The  first  one  involves  opening-dominated  inter-layer  cracking  or  “delam¬ 
ination  buckling”  (Gioia  and  Ortiz  [6];  Kadomateas  [7]).  The  second  one  involves 
shear-dominated  inter-layer  cracks  or  “shear  delaminations,”  and  often  occurs  in 
layered  materials  subjected  to  out-of-plane  impact  (Choi  et  al.  [8];  Lambros  and 
Rosakis  [9]). 

The  second  major  category  is  referred  to  as  intra-layer  failure.  There  are  three 
possible  intra-layer  failure  modes  depending  on  the  material  constitution.  The  first 


1.  Opening  dominated 
inter-layer  crack 

2.  Shear  dominated 
inter-layer  crack 


3.  Intra-layer  crack 
(matrix  crack) 


4.  Fiber  breakage 


5. fiber/matrix 
debonding 


Fig.  1  Failure  modes  for  layered  materials  based  on  material  constitution 
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Fig.  2  A  series  of  back  lit  photos  showing  the  dy  namic  failure  process  for  a  fiberglass  face  plates/ 
PVC  foam  core  sandwich  structure 


one  is  called  intra-layer  cracking  or  matrix  cracking  as  seen  in  Fig.  2.  This  type  of 
cracking  often  occurs  inside  the  matrix  of  fiber-reinforced  composites  or  within  the 
soft  core  of  sandwich  structures  [10].  It  is  also  found  in  the  form  of  tunneling  cracks 
in  thin  film/substrate  structures  [1].  Another  possible  intra-layer  failure  is  the  failure 
of  reinforcements  such  as  fiber  breakage  and  fiber  kinking  within  a  layer  (Oguni 
et  al.  [111).  The  third  intra-layer  failure  mode  is  interfacial  debonding  between  the 
matrix  material  and  the  reinforcement  such  as  debonding  between  particle/fiber  and 
matrix  occurring  within  a  constituent  layer  [12].  Indeed,  the  sequence,  nature  and 
interaction  of  these  dynamic  failure  modes  were  never  properly  clarified. 

In  perhaps  the  first  attempt  to  visualize  impact  failure  in  sandwich  structures 
used  in  Naval  applications,  Semenski  and  Rosakis  [13]  tested  thin  sections  of  such 
plate  structures  composed  of  PVC  foam  cores,  sandwiched  between  E-glass  face¬ 
plates.  A  pulsed  laser  was  used  to  illuminate  the  specimens  from  the  back  side  and 
a  high-speed  camera  recorded  the  deformation  and  failure  events.  A  sequence  of 
photographs  corresponding  to  this  process  is  shown  in  Fig.  2  together  with  the  post¬ 
mortem  picture  of  the  recorded  specimen.  As  evident  from  the  post-mortem  picture, 
there  are,  at  least,  two  types  of  failure  present.  Inter-layer  failure  demonstrates  itself 
in  the  form  of  delamination  between  the  face  plates  and  the  foam  core  at  the  vicin¬ 
ity  of  the  impact  site  and  free  edges.  On  the  opposite  side,  delamination  is  evident 
only  on  the  top  and  on  the  bottom  part  of  the  specimen.  Intra-layer  failure  in  the 
form  of  mode  I,  opening  cracks  in  the  soft  core  is  also  observed  forming  a  highly 
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symmetric  pattern.  Because  the  core  is  opaque,  the  high-speed  pictures  shown  be¬ 
low  are  of  limited  use.  What  they  show,  however,  is  the  emergence  and  propagation 
of  the  opening  intra-layer  (matrix)  cracks  inside  the  foam  core.  Indeed,  these  cracks 
seem  to  originate  at  the  fiber  glass/PVC  interface  opposite  to  the  side  of  impact, 
and  to  symmetrically  propagate  towards  the  impact  point.  These  cracks  originate  at 
the  same  location  where  the  fiberglass/PVC  delamination  terminates.  However,  the 
time  sequence  and  interaction  between  such  inter-layer  delamination  and  the  visible 
intra-layer,  opening  cracks  in  the  core  are  not  obvious.  Indeed,  the  backlit  real-time 
photographs  do  not  show  any  evidence  of  interfacial  delamination  within  the  time 
window  of  observation.  As  we  will  show  later,  this  observation  is  misleading  and  is 
due  to  the  fact  that  inter-layer  fractures  are  typically  shear-dominated.  As  such  they 
do  not  allow  for  light  to  go  through  during  the  recording  event  because  the  shear 
crack  faces  remain  in  contact  at  the  early  stages  of  this  process. 

The  inability  of  back-lit  photography  to  visualize  the  failure  process  completely, 
motivates  the  use  of  partially  transparent  model  sandwich  systems  which  allow  the 
use  of  full  field  optical  techniques  capable  of  capturing  the  nucleation  and  growth 
of  both  opening  and  shear-dominated  cracks  and  their  transition  from  one  mode  to 
the  other  [14].  For  many  complex  engineering  problems,  model  experiments  may 
prove  extremely  useful  as  intermediate  steps,  which  reveal  the  basic  physics  of 
the  problem  and  provide  relatively  straightforward  explanations  of  the  failure  pat¬ 
terns  observed  in  post-mortem  observations.  A  successful  approach  was  adopted  by 
Walter  and  Ravichandran  [15],  who  designed  a  model  aluminum/PMMA/aluminum 
specimen  to  simulate  and  visualize  the  static  debonding  and  matrix  cracking  pro¬ 
cess  in  ceramic  matrix  composites.  In  our  experiments,  we  also  adopt  the  same  idea 
and  introduce  an  appropriate  intermediate  model  configuration.  In  order  to  simulate 
the  difficult  three-dimensional  problem  of  the  out-of-plane  impact  of  real  sandwich 
structures  and  to  simultaneously  preserve  the  essence  of  the  failure  phenomena  in¬ 
volved,  we  introduce  a  two-dimensional,  plane  stress  specimen,  which  represents  a 
cross-sectional  cut  of  the  layered  material  as  illustrated  in  Fig.  3. 


Fig.  3  Model  layered 
specimens  are  idealized  cross 
sections  of  real  structure 
subjected  to  out-of-plane 
impact 
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For  this  type  of  model  specimen,  the  failure  process  is  easy  to  record,  visualize 
and  analyze.  It  is  noted  that  although  the  exact  impact  mechanics  involved  in  two 
configurations  is  not  identical  (the  real  case  is  three-dimensional  while  the  model 
specimen  is  closer  to  a  plane  stress  state),  the  mechanisms  of  stress  wave  propaga¬ 
tion  and  failure  progression  of  the  real  and  the  model  layered  materials  are  quite 
analogous.  In  designing  these  model  two-dimensional  sandwich  specimens,  it  is 
important  to  select  model  materials  whose  elastic  mismatch  is  similar  to  that  of  ma¬ 
terials  used  in  real  engineering  applications  (in  our  case  P VC/composites).  Selecting 
similar  Dundurs’  parameters  [1]  may  ensure  similarity  of  the  elasto-static  response 
for  the  interfacial  mechanics  problem.  Meanwhile,  selecting  model  material  combi¬ 
nations  with  similar  ratios  of  wave  speeds  of  two  constitution  materials  to  the  real 
structure  is  perhaps  the  most  important  consideration  in  the  dynamic  case,  where 
timing  of  events  and  stress  intensity  are  governed  by  the  constituent  material  wave 
speeds.  Also,  the  ratio  of  inter-layer  and  intra-layer  strengths  (or  fracture  tough¬ 
nesses)  is  important.  These  three  issues  provide  sets  of  similarity  rules  to  connect 
the  real  structures  to  our  models  experiments.  As  schematically  shown  in  Fig.  3, 
matrix  cracks  and  delamination  are  the  two  major  impact  failure  modes  in  sandwich 
structures  and  layered  composites.  At  some  intersection  points,  matrix  cracks  and 
delamination  are  connected  as  also  seen  in  the  post-impact  picture  of  Fig.  2.  One 
frequently  asked  question  in  the  literature  is  whether  the  matrix  cracks  lead  to  the 
delamination  or  the  delamination  happens  first  and  subsequently  kinks  into  the  ad¬ 
jacent  layer  inducing  the  matrix  crack.  This  is  a  typical  problem  of  sequence  and 
failure  mode  transition  identification. 

Since  the  nature  and  origin  of  such  failure  mechanisms  can  only  be  theorized 
by  post-mortem  observations,  the  necessity  of  full-field  real-time,  high-speed  mea¬ 
surements  becomes  obvious.  To  this  effect,  the  objectives  of  the  current  work  are 
to  conduct  systematic  experimental  studies  of  the  time  evolution  and  nature  of  dif¬ 
ferent  dynamic  failure  modes  and  to  investigate  their  interactions.  Through  these 
model  experiments,  we  try  to  identify  the  basic  physical  phenomena,  and  to  provide 
guidance  for  theoretical  models  and  much  needed,  real-time,  validation  of  numer¬ 
ical  codes.  To  make  this  comparison  more  meaningful,  we  choose  model  material 
combinations  that  have  the  ratios  of  wave  speeds  very  close  to  those  used  in  real 
sandwich  structures. 


2  Experimental  Procedure 
2.1  Materials  and  Specimens 

Two  kinds  of  materials  were  used  in  the  experiments  described  below.  A  4340- 
carbon  steel  was  employed  to  simulate  the  stiff  and  strong  fiberglass  faceplates  of 
sandwich  structures.  The  polymeric  material,  which  was  used  to  simulate  the  weak 
core  layer,  such  as  the  PVC  foam  core  or  balsa  wood  in  sandwich  structures  or 
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the  90°  plies  in  cross-ply  laminates,  is  Homalite-100.  The  adhesive  used  to  bond 
the  metal/polymer  interface  is  Weldon- 10.  The  detailed  properties  of  this  adhesive 
and  the  effect  of  interfacial  strength  variation  on  dynamic  failure  mode  selection 
are  reported  by  Xu  and  Rosakis  [16,  17].  The  shear  wave  speed  is  an  important 
parameter  in  this  investigation.  The  shear  wave  speed  ratio  for  the  core  and  faceplate 
is  3.2  for  typical  E-glass/PVC  sandwich  structures  of  the  type  that  have  recently 
been  used  in  construction  of  full-scale  composite  ships.  Details  on  the  complete  set 
of  physical  and  constitutive  properties  for  E-glass  composite  materials  have  recently 
been  discussed  by  Oguni  et  al.  [11],  For  comparison,  the  same  shear  wave  speed 
ratio,  for  the  idealized  steel/Homalite  model  sandwiches,  is  about  2.7  based  on  the 
data  from  Ref.  [14],  Although  the  absolute  values  of  these  constituent  properties 
are  very  different  in  the  “idealized”  versus  the  “real”  solids,  the  idealized  material 
combinations  have  been  chosen  in  such  a  way  as  to  have  a  shear  wave  speed  ratio 
that  is  very  similar  to  its  real  sandwich  counterparts. 

As  shown  in  Fig.  4,  three  different  types  of  model  sandwich  specimen  geometries 
were  designed  and  tested.  Type  A  specimens  have  equal  layer  widths  and  involve 
two  different  materials.  They  contain  two  metal  layers  with  one  polymer  layer  sand¬ 
wiched  between  them.  Type  B  specimens  involve  two  thin  metal  layers  (faceplates) 
and  one  polymer  layer.  This  type  of  specimens  is  quite  similar  in  geometry  (ratio  of 
core  to  face  plate  thickness)  to  realistic  sandwich  plates  used  in  engineering  appli¬ 
cations.  The  only  difference  between  type  C  and  type  B  specimens  is  their  lengths. 
Type  C  specimens  are  twice  as  long  as  type  B  specimens.  The  purpose  of  type  C 
specimens  is  to  explore  the  impact  failure  patterns  with  least  edge  effect  present  in 
the  time  scale  of  the  failure  process.  All  three  types  of  specimens  have  the  same 
out-of-plane  thickness  of  6.35  mm. 
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Fig.  4  Model  specimens  simulating  sandwich  structures  ( shaded  layers  -  metals;  transparent 
layers  -  polymers) 
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The  majority  of  experiments  in  this  investigation  were  performed  using  dynamic 
photoelasticity.  This  classical  method  has  recently  found  a  lot  of  new  applications 
such  as  study  of  the  dynamic  fracture  processes  in  functionally  gradient  materials 
(FGMs)  [18].  The  Coherent  Gradient  Sensing  (CGS)  method  [19]  was  also  used  in  a 
small  number  of  cases.  A  schematic  of  the  dynamic  photoelasticity  setup  used  here 
is  given  in  Fig.  5.  Two  circular  polarizer  sheets  were  placed  on  either  side  of  the 
specimen.  A  coherent,  monochromatic,  plane  polarized  laser  light  output  was  col¬ 
limated  to  a  circular  beam  of  100  mm  in  diameter.  The  laser  beam  was  transmitted 
through  the  specimen  and  the  resulting  fringe  pattern  was  recorded  by  the  high¬ 
speed  camera.  A  Cordin  model  330A  rotating-mirror  type  high-speed  film  camera 
was  used  to  record  the  images.  During  the  impact  test,  a  projectile  was  fired  by  the 
gas  gun  and  impacted  the  specimen  center.  The  generation  of  isochromatic  fringe 
patterns  is  governed  by  the  stress-optic  law.  For  the  case  of  monochromatic  light, 
the  condition  for  the  formation  of  fringes  can  be  expressed  as  [20]: 


(Ti  —  ct2  = 


Nfn 

h 


(1) 


where  d\  —  <r2  is  the  principal  stress  difference  of  the  thickness  averaged  stress 
tensor.  fa  is  the  material  fringe  value,  N  is  the  isochromatic  fringe  order  and  h  is 
the  specimen  thickness.  The  isochromatic  fringe  patterns  observed  are  proportional 
to  contours  of  constant  maximum  in-plane  shear  stress,  rmax  =  (cti  —  <t2)/2. 


Fig.  5  Schematic  of  the  dynamic  photoelasticity  and  high-speed  photography  setup 
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3  Results  and  Discussion 

3.1  Failure  Process  of  Short  Model  Sandwich  Specimens 
with  Equal  Layer  Widths 

The  diameter  of  the  laser  beam  was  100  mm;  however,  the  maximum  length  of  the 
zone  that  had  to  be  investigated  was  254  mm  long.  In  order  to  observe  all  possible 
dynamic  failure  modes  present  in  each  case,  the  field  of  view  had  to  be  moved 
from  one  location  to  another  for  each  specimen  configuration  under  the  same  impact 
condition.  Figure  6  presents  a  series  of  photoelastic  images  of  the  Homalite  core 
layer  of  a  type  A  specimen. 

The  dark  circular  spot  at  the  upper  right  comer  is  a  scaling  mark  (diameter 
6.35  mm)  bonded  on  the  specimen.  The  thin  horizontal  dark  line,  seen  around  the 
center  of  every  image,  is  the  streak  line  of  the  camera.  This  line  provides  a  station¬ 
ary  reference  when  the  whole  specimen  moves  during  the  impact  process.  At  first, 
the  field  of  view  was  centered  on  the  middle  of  the  specimen  because  it  was  close  to 
the  impact  position  and  failure  was  expected  to  initiate  from  this  zone.  As  shown  in 
Fig.  6b,  about  158  (is  after  impact,  two  inter-layer  cracks  at  the  lower  interface  en¬ 
tered  the  field  of  view  and  propagated  towards  the  specimen  center.  Later  on  (around 
1 82  |xs),  the  two  inter-layer  cracks,  identified  by  the  moving  concentration  of  fringes 
at  their  tips,  are  seen  to  meet  each  other  in  Fig.  6d.  Similar  to  shear-dominated  inter¬ 
facial  cracks  in  bimaterials  [21],  those  inter-layer  cracks  are  also  shear-dominated. 
Because  the  Homalite  and  steel  layers  are  still  in  contact  up  to  that  time,  no  visual 
evidence  of  decohesion  is  apparent  in  the  images,  although  these  cracks  have  al¬ 
ready  broken  the  interface  in  a  combination  of  compression  and  shear.  After  these 
two  inter-layer  cracks  meet  at  the  center,  a  bright  gap  between  the  Homalite  and 
steel  layers  can  be  seen  to  appear  in  Fig.  6e.  Along  this  clearly  opened  interface 
and  on  the  Homalite  side,  two  Rayleigh  surface  waves  are  now  seen  to  propagate, 
originating  from  the  center  and  moving  outwards,  towards  the  specimen  edges. 

Perhaps  the  most  interesting  conclusion  deriving  from  the  sequence  shown  in 
Fig.  6  is  the  fact  that  delamination  did  not  initiate  in  the  interface  directly  above  the 
impact  point  but  did  so  outside  our  central  field  of  view  at  two  symmetric,  off-axis, 
locations  along  the  lower  interface.  In  order  to  discover  the  location  of  crack  nu- 
cleation,  we  must  move  our  field  of  view  off  the  specimen  center  to  investigate  the 
origins  of  these  inter-layer  cracks.  In  order  to  conclusively  identify  the  origins  of 
the  upper  and  lower  inter-layer  cracks,  the  field  of  view  was  once  more  moved  to  the 
specimen  edge  as  shown  in  Fig.  7a.  After  impact  at  the  specimen  center,  the  stress 
waves  in  the  bottom  steel  layer  propagated  towards  the  edge  creating  a  visible  head 
wave  structure  on  the  lower  wave  speed  polymer  side  (see  Fig.  7b).  Right  after  the 
stress  wave  reached  the  free  edge,  due  to  the  existence  of  a  stress  singularity  at  the 
bimaterial  comer  [22,23],  an  inter-layer  crack  initiated  at  the  lower  interface  as  seen 
in  Fig.  7c.  This  crack  propagated  towards  the  specimen  center.  After  around  160  (is, 
another  inter-layer  crack  initiated  at  the  upper  interface  also  moving  towards  the 
center.  This  upper  inter-layer  crack  soon  kinked  into  the  core  layer  and  branched 
into  a  fan  of  multiple  mode-I  intra-layer  cracks.  Crack  branching  as  reported  by 
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Fig.  6  Early  stages  of  the  failure  process  of  a  three-layer  specimen  with  equal  layer  widths.  The 
central  field  of  view  reveals  the  early  occurrence  of  shear-dominated  delamination  at  the  lower 
interface.  At  later  times  (Fig.  fie),  the  debonding  becomes  clearly  visible  at  the  lower  interface 


previous  researchers  often  initiates  when  a  crack  in  a  homogeneous  solid  reaches 
high  fractions  of  the  shear  wave  speed  [24],  for  example,  30-40%  shear  wave  speed 
of  Homalite-100. 

Based  on  experimental  observations  from  different  fields  of  view,  the  major 
dynamic  failure  modes  and  sequence  in  model  three-layer  materials  can  be  sum¬ 
marized  in  Fig.  8.  After  the  stress  wave  reaches  the  free  edges,  two  shear-dominated 
inter-layer  cracks  initiate  and  propagate  towards  the  specimen  center.  These  shear 
cracks  separate  the  whole  lower  interface  and  a  Rayleigh  surface  wave  forms  on  the 
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Fig.  7  Edge  view  of  damage  evolution:  delaminations  are  shown  to  form  at  the  intersection  of, 
first  lower,  and  then  upper  interfaces  with  the  specimen  edge.  A  fan  of  kinked  intra-layer  cracks 
branches  into  the  core  layer  from  the  upper  interface 


separated  free  surface.  At  a  later  stage,  inter-layer  cracks  also  originate  from  the  up¬ 
per  interface  at  the  free  edge  and  travel  towards  the  specimen  center.  However,  these 
upper  inter-layer  cracks  soon  kink  into  the  core  layer  to  form  opening-dominated 
intra-layer  cracks.  Under  certain  circumstances  (e.g.,  if  the  core  material  is  very 
brittle),  such  kinked  cracks  may  also  branch  into  a  fan  of  multiple  branches  frag¬ 
menting  the  core.  The  model  experiments  described  here  seem  to  capture  the  basic 
nature  of  the  post-mortem  impact  failure  modes  observed  in  real  sandwich  struc¬ 
tures.  Indeed,  the  kinked  matrix  crack  of  the  core  layer  of  the  glass  fiber/foam 
core  sandwich  shown  in  Fig.  8e  seems  to  follow  the  same  initiation  and  propagation 
process  as  the  kinked  intra-layer  crack  in  the  model  three-layer  specimens  schemat¬ 
ically  shown  in  Fig.  8d. 


3.2  Failure  Process  in  Long  Model  Sandwich  Specimens 

In  type  A  specimens,  inter-layer  cracks  always  initiated  from  the  specimen  free 
edges  due  to  the  bi-material  stress  concentration  at  such  locations.  In  order  to  study 
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3  Lower  interfacial  debonding  b  Completely  debonding 


V  1 


Rayleigh  wave  traveling  on  the  opening  interface 


Fig.  8  Conceptual  summary  of  typical  failure  modes  and  sequence  (a-d)  in  a  short  three-layer 
specimen  with  equal  layer  widths,  (e)  Is  a  comparison  with  a  post-impact  picture  of  a  real  sandwich 
specimen 

the  impact  damage  modes  and  failure  sequence  in  either  very  large  structures  or  ones 
that  are  clamped  along  the  edges,  our  model  specimens  featured  long  specimens. 
These  specimens  were  long  enough  such  that  any  damage  from  the  edges,  such  as 
inter-layer  cracks  induced  by  the  edge  effect,  arrived  in  the  area  of  observation  long 
after  the  local  damage  sequence  had  been  completed  near  the  impact  point.  To  study 
this  effect,  we  tested  long  sandwich  style  specimens  of  type  C  (see  Fig.  9). 

As  shown  in  Fig.  9d,  at  79  p  s  after  impact,  an  inter-layer  crack  tip  is  seen  at  the 
lower  interface  (fringe  concentration  within  the  dashed  circle).  This  crack  is  similar 
in  nature  to  our  previously  observed  inter-layer  fractures  but  has  not  originated  at 
the  specimen  free  edge,  which  for  type  C  specimens  is  far  away  from  our  field  of 
view.  Indeed,  if  this  crack  originated  from  the  specimen  free  edge,  it  would  take  at 
least  150  (is  to  enter  the  field  of  view.  Closer  scrutiny  reveals  that  this  crack  origi¬ 
nates  from  a  much  closer  location  to  the  impact  point.  This  location  is  marked  here 
by  the  circle  in  Fig.  9c  within  which  a  concentration  of  photoelastic  fringes  points 
to  the  concentration  of  shear  stresses  that  is  responsible  for  its  nucleation.  Indeed 
the  crack  nucleates  at  a  location  where  the  inter-layer  shear  stress  reaches  a  local 
maximum,  whose  value  equals  the  shear  strength  of  the  bond.  To  rationalize  this, 
one  should  consider  the  symmetry  of  our  impact  configuration  and  recall  the  strong 
wave  speed  mismatch  between  the  lower  faceplate  and  the  core  material.  The  shear 
stress  component  on  at  the  specimen  centerline  will  always  vanish  because  of  this 
symmetry  but  is  expected  to  anti-symmetrically  increase  away  from  the  centerline 
as  compressive  waves  begin  to  spread  along  the  steel  faceplate.  As  a  result  of  these 
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Fig.  9  Nucleation  of  an  intersonic  inter-layer  crack  at  the  vicinity  of  the  impact  area  of  a  long 
model  sandwich  specimen 


observations,  the  scenario  that  seems  to  be  emerging  is  as  follows:  shear-dominated 
cracks  are  generated  at  two  points  to  the  right  and  left  of  the  center  line  and  more 
backwards  towards  the  impact  point.  A  series  of  photographs  confirming  the  exis¬ 
tence  of  an  inter-layer  crack  coming  from  the  right-hand  side  of  the  impact  point 
is  shown  in  Fig.  9.  Indeed  Fig.  9c  corresponds  to  the  nucleation  of  this  crack  while 
Fig.  9d-f  confirm  its  high-speed  motion  towards  the  impact  site.  As  this  inter-layer 
crack  and  its  symmetric  companion  from  the  left  meet  above  the  impact  point,  they 
create  a  central  shear  delamination  between  the  core  and  the  bottom  faceplate.  The 
speed  of  this  crack  is  very  high  as  evident  from  the  shear  shock  wave  that  appears 
as  a  dark  inclined  line  radiating  from  its  moving  tip  (Fig.  9e,  f)- 

Figure  10  corresponds  to  another  impact  experiment  which  featured  the  same 
load  condition  as  shown  in  Fig.  9.  The  end  point  of  the  central  delamination 
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Fig.  10  Local  view  of  the  post-mortem  damage  in  a  long  sandwich  specimen  (a)  and  high-speed 
snap  shot  capturing  the  formation  of  intra-layer  cracks  (b) 


described  above  is  denoted  by  A.  Yu  et  al.  [25]  show  that  the  inter-layer  normal 
stresses  near  the  impact  point  are  compressive  while  the  inter-layer  shear  stresses 
exhibit  two  clear  peak  values  (same  magnitude)  moving  symmetrically  away  from 
the  centerline.  Point  A  in  Fig.  10  is  corresponding  to  the  old  location  of  the  maxi¬ 
mum  inter-layer  shear  stress.  This  point  can  now  acts  as  a  stress  concentration  from 
which  further  damage  (to  the  core  as  well  as  to  the  rest  of  the  interface)  will  subse¬ 
quently  evolve.  Indeed  as  seen  in  Fig.  10,  intra-layer  cracks  now  are  generated  and 
propagated  into  the  core  (along  AC),  also  accompanied  by  a  new  inter-layer  debond 
(along  AB)  also  originating  at  point  A.  The  high-speed  snapshot  that  appears  in  the 
same  photograph  confirms  this  scenario. 

Figure  1 1  summarizes  the  proposed  failure  evolution  sequence  for  the  long  sand¬ 
wich  style  specimens  described  above.  One  point  that  should  be  made  clear  here  is 
that  following  the  formation  of  the  central  (shear)  delamination,  the  choice  of  the 
inclination  angle  f)  and  the  possibility  of  further  delamination  along  the  bottom  inter¬ 
face  depend  on  the  impact  speed  and  on  the  relative  values  of  the  matrix  material  and 
interfacial  bond  strengths.  The  same  is  true  for  the  exact  locations  of  points  A  and  B. 
However,  we  expect  that  if  impact  speeds  are  high  enough  to  promote  this  localized 
failure  mode,  the  general  features  described  here  will  continue  appearing  even  as 
the  projectile  speed  increases  further.  For  the  initiation  of  intra-layer  cracks  (matrix 
cracks),  previous  researchers  theorized  that  such  cracks  initiated  from  the  center 
of  the  weak  layer  and  propagated  toward  adjacent  interfaces  to  lead  to  inter-layer 
cracks  or  delaminations  [26].  However,  no  real-time  experimental  evidence  was 
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Fig.  11  Failure  sequence  observed  in  long  sandwich-style  specimens  with  minimal  edge  effects 


ever  observed  to  support  such  a  scenario.  Here,  we  clearly  show  that  the  intra-layer 
cracks  always  initiate  at  the  interfaces  following  the  shear-dominated  delamination, 
which  kinks  into  the  core  layer  resulting  in  intra-layer  cracking. 


3.3  Effect  of  Impact  Speeds 

In  Section  3.1,  we  investigated  the  nature  and  sequence  of  failure  mechanisms  in 
relation  to  model  sandwich  specimens  subjected  to  an  impact  speed  of  33  m/s,  as 
seen  in  Fig.  6.  This  impact  speed  situation  will  be  taken  as  the  baseline  for  our  com¬ 
parisons.  Figure  12  describes  an  experiment  of  the  same  geometry  that  corresponds 
to  an  impact  speed  is  45  m/s  [16].  As  seen  from  Fig.  12b,  two  inter-layer  cracks 
appeared  at  the  lower  interface  and  propagated  towards  the  center,  racing  towards 
each  other  with  intersonic  speeds.  At  a  later  time,  inter-layer  cracks  at  the  upper 
interface  also  appeared  propagating  towards  the  center  (Fig.  12c).  The  locations  of 
these  four  inter-layer  cracks  (two  at  the  top  and  two  at  the  lower  interfaces)  are  in¬ 
dicated  by  the  white  arrows.  As  clearly  seen  from  Fig.  12e,  intra-layer  damage  also 
spreads  from  the  interface  in  to  the  core  in  the  form  of  a  periodic  series  of  mode-I 
cracks  inclined  at  a  small  angle  to  the  vertical  axis.  These  cracks  are  nucleated  at 
the  upper  interface  at  locations  that  are  behind  the  horizontally  moving  inter-layer 
shear  crack.  Their  nucleation  and  growth  result  in  the  eventual  fragmentation  of  the 
specimen  core.  The  inter-layer  cracks  propagating  at  the  lower  bimaterial  interface 
and  facing  towards  each  other  in  Fig.  12d-e  feature  clearly  formed  shock-like  or 
Mach-like  discontinuities  (shear  shock  waves)  which  are  emitted  from  their  crack 
tips.  These  discontinuities  in  photoelastic  patterns  represent  traveling  discontinuities 
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Fig.  12  Growth  of  four  inter-layer  cracks  at  the  center  of  a  three-layer  specimen 

in  maximum  shear  stress  and  are  clear  proofs  that  crack  tips  have  exceeded  the 
shear  wave  speed  of  Homalite  [19].  These  shock  waves  formed  a  clear  testimony  to 
the  intersonic  nature  of  the  inter-layer  crack  growth  even  before  any  detailed  crack 
measurement  was  ever  attempted.  The  crack  speed  history  for  the  lower,  right  inter¬ 
layer  crack  is  plotted  in  Fig.  13  as  a  function  of  distance  from  the  free  edge.  The 
figure  shows  that  the  crack  speed  of  the  higher  impact  speed  case  (45  m/s)  is  always 
higher  than  the  baseline  equivalent  remaining  always  intersonic  within  the  window 
of  observation. 


3.4  Dynamic  Failure  Mode  Transition 

From  the  above  experimental  observation,  we  find  that  different  dynamic  failure 
modes  and  their  transition  are  very  interesting.  As  seen  in  Fig.  lib,  an  interface 
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Fig.  13  Comparison  of  crack  speed  distributions  of  two  identical  specimens  subjected  to  different 
impact  speeds.  The  dash  line  is  the  dynamic  shear  wave  speed  of  the  Homalite-100 
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Fig.  14  Common  failure  modes  when  a  crack  encounters  an  interface 


crack  kinks  from  the  interface  to  form  a  matrix  or  intra-layer  crack.  While  this 
matrix/intra-layer  will  deflect  at  the  second  interface  to  form  a  inter-layer  crack 
again.  Generally,  when  a  crack  propagates  in  elastic  solids  and  encounters  an  in¬ 
terface,  one  of  the  three  situations  may  occur  as  seen  in  Fig.  14:  (a)  after  the  crack 
reaches  the  interface,  it  kinks  out  of  its  original  path  and  continues  to  propagate 
along  the  interface  [27].  This  phenomenon  is  often  called  “crack  kinking  or  deflec¬ 
tion”  [28];  (b)  the  crack  penetrates  the  interface  and  continues  to  propagate  along 
its  original  path,  i.e.,  crack  penetration  [29,30];  (c)  early  interface  debonding  initi¬ 
ates  before  the  incident  crack  reaches  the  interface,  or  it  refers  to  the  “Cook-Gordon 
mechanism”  [31-33].  In  the  open  literature,  efforts  have  been  primarily  focused  on 
analyzing  the  first  two  cases,  crack  kinking  and  crack  penetration  [1],  The  energy 
release  rate  ratios  of  the  incident  and  kinked  interfacial  cracks,  and  the  fracture 
toughness  ratios  of  the  matrix  material  and  the  interface  are  identified  as  major 
parameters  to  govern  crack  deflection/penetration  [34], 

Recently.  Xu  et  al.  [30]  experimentally  and  analytically  studied  the  dynamic 
crack  deflection/penetration  phenomena.  They  also  presented  an  energy-based 
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criterion  to  investigate  the  competition  between  the  dynamic  crack  penetration  and 
deflection.  The  criterion  for  interfacial  dynamic  crack  kinking  is  a  ratio  comparison 
between  two  energy  release  rates  for  the  deflected/kinked  crack  along  the  interface, 
and  the  incident  crack  inside  the  matrix  (driving  force)  and  the  fracture  toughness 
values  (material  resistance): 


Gd0.Vl)  Tf  (V;) 

GfM  ~  r, '}'•(»,) 

where  is  the  crack  kinking  angle  or  interface  angle,  and  vj  and  v’2  are  the  incident 
parent  crack  speed  and  the  kinked  daughter  crack  speed  respectively  (Freund  [401). 
r">0  is  the  dynamic  fracture  toughness  of  the  matrix  material  and  r^r(v2)  is  the 
interfacial  fracture  toughness.  If  the  left  side  is  less  than  the  right  side,  a  dynamic 
crack  penetration  will  occur. 

However,  in  order  to  apply  the  energy-based  criterion,  putative  crack  deflection 
and  crack  penetration  lengths  are  needed  to  evaluate  the  corresponding  energy  re¬ 
lease  rates.  Several  researchers  [29,  35]  have  demonstrated  that  the  two  putative 
lengths  have  a  significant  effect  on  the  energy  release  rate  ratios,  and  sometimes  the 
energy-based  criterion  fails  to  predict  the  crack  deflection  or  crack  penetration.  For 
these  cases,  “Cook-Gordon  mechanism"  provides  an  alternative  explanation  since  a 
crack  may  not  kink  right  after  it  reaches  the  interface  as  shown  in  Fig.  14a.  The  case 
(interface  debonding  before  kinking)  shown  in  Fig.  14c  is  quite  possible. 


3.5  Dynamic  Interface  Debonding  Ahead 

of  a  Main  Incident  Crack 

In  Fig.  14c,  correlations  of  the  fracture  mechanics  parameters  of  the  kinked  inter¬ 
facial  crack  and  the  incident  crack  are  not  easy  to  obtain.  Therefore,  in  order  to 
model  the  “Cook-Gordon  mechanism",  we  tend  to  use  strength-based  criteria  to  pre¬ 
dict  interfacial  debonding  initiation  only  (rather  than  crack  growth)  induced  by  an 
incident  crack.  In  terms  of  the  dynamic  “Cook-Gordon  mechanism",  only  Needle- 
man  and  co-workers  [36-38]  have  simulated  this  problem  using  a  cohesive  element 
model.  In  their  model,  an  artificial  initial  flaw  was  introduced  so  they  assumed  some 
material  properties  for  predictions.  In  our  investigation,  a  strength  criterion  with  di¬ 
rect  interfacial  strength  measurements  will  be  used  to  predict  the  critical  distance 
rc  of  the  incident  crack  tip  to  the  intersection  point  of  the  incident  crack  path  and 
the  interface.  Indeed,  our  work  will  be  complementary  to  Needleman’s  work,  since 
our  work  aims  to  predict  interfacial  debonding  initiation,  while  their  efforts  were 
focused  on  simulating  the  late  interfacial  crack  propagation  after  crack  initiation. 

Figure  15  shows  a  series  of  photoelasticity  snap  shots  following  impact  of  a 
bonded  Homalite  specimen  [33].  In  all  experiments,  the  projectile  impacted  the  cen¬ 
ter  of  the  bottom  layer  on  a  steel  buffer  as  shown  in  Fig.  15a.  Figure  15b  shows  a 
fan  of  mode  I  cracks  (symmetric  fringe  patterns)  appearing  from  the  upper  free  edge 
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Fig.  15  Dynamic  failure  process  in  a  two-layer  specimen  showing  the  interaction  of  a  fan  of  mode 
1  incident  cracks  and  the  resulting  interfacial  crack.  The  thin  horizontal  line  is  the  weak  interface 


at  approximately  93.8  |is  after  impact.  After  impact,  the  longitudinal  compressive 
stress  wave  traveled  from  the  lower  impact  side  towards  the  upper  free  edge. 

This  compressive  stress  w'ave  reflected  from  this  edge  as  a  tensile  w;ave  and  its 
intensity  was  sufficient  to  nucleate  a  fan  of  branched  cracks  from  the  free  edge.  As 
time  goes  on  (Fig.  15c— f),  the  nucleated  fan  of  cracks  widens  significantly  by  pro¬ 
ducing  a  multiplicity  of  both  successful  and  unsuccessful  branches,  some  of  which 
move  tow'ards  the  still  coherent  interface  (for  a  discussion  of  crack  branching  phe¬ 
nomenon  in  bulk  Homalite,  see  Ravi-Chandar  and  Knauss  [24],  The  average  speed 
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of  these  locally  mode  I,  branched  cracks  is  0.41  Cs,  which  is  the  branching  speed 
in  bulk  Homalite.  Well  before  the  branched  cracks  reached  the  interface,  a  central 
inter-layer  crack  was  nucleated  at  the  intersection  of  the  specimen  center  line  and  the 
bond  line  as  seen  in  Fig.  15c.  This  interfacial  crack  propagated  in  both  directions  off 
the  center  as  shown  in  Fig.  15d.  At  the  specimen  centerline,  the  shear  stresses  vanish 
because  of  symmetry.  As  a  result,  the  nucleated  inter-layer  crack  is  initially  and,  for 
a  very  short  time,  mode  I  dominated.  Its  nucleation  is  induced  by  the  stress  field 
produced  by  the  fan  of  branched  cracks  approaching  the  interface.  As  this  crack 
spreads  symmetrically,  opening  up  the  interface  (see  distinct  evidence  of  decohe¬ 
sion  in  Fig.  15f,  the  fan  of  branched  cracks  decelerates  and  arrests  just  before  these 
cracks  reach  the  decohered  interface.  The  above  described  scenario  is  perhaps  the 
first  real-time  visualization  of  the  dynamic  equivalent  of  the  “Cook-Gordon  Mech¬ 
anism”  describing  the  remote  decohesion  of  an  interface  due  to  the  approach  of 
a  matrix  (intra-layer)  crack.  As  the  interfacial  crack  spreads  away  from  the  speci¬ 
men  centerline,  it  almost  immediately  encounters  increasing  amount  of  interfacial 
shear  stress,  which  quickly  converts  it  to  a  mixed-mode  and  eventually  to  a  mode  II 
dominated  crack.  Unlike  propagating  cracks  in  bulk  Homalite,  interfacial  cracks  are 
constrained  to  propagate  along  the  weak  interface  and,  as  a  result,  they  can  do  so 
under  mixed-mode  or  primarily  mode  II  conditions.  They  can  also  propagate  at  very 
high  (even  intersonic)  speeds  compared  to  their  bulk  (intra-layer)  counterparts.  In 
order  to  predict  the  dynamic  crack  initiation,  Wang  and  Xu  [27]  proposed  a  strength 
criterion 

(02J0,)1  +  (a[2/as)2  =  1  (3) 

where  a,,as  are  the  interface  tensile  and  shear  strengths,  a'n , a[2. a22  are  stresses 
acting  on  the  interface  which  are  induced  by  an  incident  dynamic  crack  close  to  the 
interface.  The  stress  field  of  a  steady  mode  I  crack  inside  a  linear  elastic  solid  is 
given  by  a  well-known  form  [39,40]: 

°\i  =  +  TSnSjt  +  0(1)  (4) 

■Jlnr 

where  K/(t)  is  the  dynamic  stress  intensity  factor  of  the  mode  I  crack  as  a  func¬ 
tion  of  time  t\T  is  a  non-singular  term,  which  is  called  “the  T-stress”  or  aox  [20]; 
0(1)  represents  higher  order  terms;  the  functions  2,y(#,  v)  that  represent  the  angu¬ 
lar  variation  of  stress  components  for  an  instantaneous  crack  tip  speed  v.  For  the 
kinked/deflected  interface  crack,  it  is  a  mixed-mode  crack  with  a  stress  field: 


_  _  /  .  _// 
aij  ~  aij  "f  °ij 


=  T(  (fl.  V)  +  TS„S„  + 

y/2 nr 


Knit) 

sflirr 


y ia(o,  v)  +  od) 


(5) 


where  Kj/(t)  is  the  dynamic  stress  intensity  factor  of  the  mode  II  crack  as  a  func¬ 
tion  of  time  t  and  the  functions  2^(0,  v)  represents  the  angular  variation  of  stress 
components  for  an  instantaneous  crack  tip  speed  v. 
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3.6  New  Progress  on  Dynamic  Crack  Branching 

As  observed  in  the  above  dynamic  failure  experiments,  we  find  that  dynamic  crack 
branching  is  an  important  dynamic  failure  mode  for  brittle  materials,  and  for  ductile 
materials  subjected  to  certain  conditions.  This  phenomenon  has  received  extensive 
attention  in  the  past  decades.  However,  previous  major  efforts  have  primarily  fo¬ 
cused  on  analytical  and  numerical  studies  [41-43].  Very  few  new  experimental 
results  were  available  to  verify  predictions  or  to  provide  guidance  for  modeling 
[24,44],  Some  important  issues,  such  as  the  crack  speed  change  before  and  after 
branching,  effect  of  dynamic  loading  rate  on  the  crack  branching,  and  crack  branch¬ 
ing  induced  by  stress  wave  loading  are  still  open.  Therefore,  we  try  to  focus  on 
how  to  realize  various  forms  of  dynamic  crack  branching  under  a  variety  of  condi¬ 
tions  [45].  The  objectives  are  to  elucidate  a  series  of  new  phenomena;  and  to  provide 
guidance  for  developing  theoretical  models  and  validating  numerical  simulations. 


3.6.1  Special  Experiments  for  Dynamic  Crack  Kinking  and  Branching 

Two  kinds  of  polymeric  materials  were  used  in  conjunction  with  two  kinds  of  optical 
diagnostic  techniques.  Homalite-100  was  chosen  for  the  photoelasticity  experiments 
while  PMMA  was  used  in  the  Coherent  Gradient  Sensing  experiments.  Various 
types  of  specimens  were  designed  and  some  of  them  had  pre-notches  with  differ¬ 
ent  radii.  One  major  specimen  used  in  this  investigation  was  a  novel  wedge-loaded 
specimen,  which  was  designed  to  produce  a  single,  straight  dynamic  crack  as  shown 
in  Fig.  16a.  An  aluminum  wedge  was  inserted  into  a  pre-notch  and  impacted  by  a 
projectile,  causing  the  wedge  to  open  the  notch  faces  thus  producing  a  single  mode 
I  crack  [30].  The  notch  tip  was  cut  using  a  diamond  afering  blade  (Buehler,  Series 
15  LC).  A  strain  gauge  was  bonded  onto  the  wedge  to  trigger  the  high-speed  camera 
and  laser  system. 


3.6.2  Dynamic  Crack  Branching  and  Kinking  from  a  Weak  Interface 

As  shown  in  Fig.  16,  the  in-plane  Homalite  specimen  dimensions  were  457  mm 
long,  254  mm  wide  and  the  plate  thickness  was  9.5  mm.  In  this  photoelasticity  ex¬ 
periment,  the  initial  notch  radius  was  0.127  mm  (0.005  in.).  For  a  low  impact  speed, 
V  =  19  m/s,  only  a  pure  mode  I  crack  initiated  from  the  notch.  As  seen  in  Fig.  16b. 
we  created  an  artificial  interface  (an  inclined  thin  line)  in  front  of  the  horizontally 
propagating  crack.  The  incident  mode  I  crack  approached  the  interface  at  about 
151  [is  after  impact  and  transitioned  into  a  mixed-mode  interfacial  crack.  A  ver¬ 
tical  line  appearing  in  every  image  is  the  camera  streak  line,  which  was  used  for 
positioning  and  reference  purposes. 

At  approximately  1 77  [is,  this  mixed-mode  interfacial  crack  kinked  into  the  right 
side  of  the  interface.  A  significant  caustics  (or  shadow  spot)  is  seen  in  Fig.  16d  to 
show  the  mode  I  nature  of  the  kinked  crack.  The  speed  of  the  kinked  crack  was 
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Fig.  16  Formation  and  propagation  of  branched  cracks  kinked  from  an  interfacial  crack 


high  enough  to  induce  multiple  branches,  which  are  visible  in  Fig.  I6e,  f.  Crack 
length  and  speed  records  of  the  main  crack  and  the  branched  cracks  are  presented  in 
Fig.  17.  Differentiation  of  the  crack  length  record  furnishes  the  tangential  crack  tip 
speed  before  and  after  crack  branching.  Since  the  differentiation  process  is  based  on 
a  three-point-fitting  procedure  of  the  crack  length  history,  the  exact  crack  speed  at 
the  crack  branching  could  not  be  obtained.  However,  it  is  interesting  to  notice  that 
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Fig.  17  Crack  length  history  and  fitted  crack  tip  speeds  of  the  main  crack  and  the  branched  cracks 


the  main  crack  tip  speed  is  almost  equal  to  the  branched  crack  tip  speed  under  the 
current  time  resolution  (2.6  ps/frame).  The  crack  branching  speed  was  about  28% 
of  the  shear  wave  speed  of  Homalite  at  the  high  strain  rate. 


3.6.3  Dynamic  Crack  Branching  Initiated  from  a  Notch  Subjected 
to  High  Impact  Loading 

If  we  kept  all  other  conditions  from  the  previous  case  but  raised  the  impact  speed 
to  30  m/s  (58%  increase),  a  mode  I  crack  was  observed  at  first  as  shown  in  Fig.  18b. 
However,  this  main  crack  soon  branched  into  two  cracks  as  seen  in  Fig.  18c.  At 
a  later  time,  one  branched  crack  (upper  branch)  generated  two  new  sub-branching 
cracks  shown  in  Fig.  18d.  Previous  experiments  reported  that  the  branching  angle 
was  less  than  45°  [20,46].  However,  in  our  experiment,  the  first  branching  angle 
is  around  45°  and  the  second  branching  angle  is  about  67°,  which  is  approximately 
the  theoretical  branching  angle  (60°)  first  reported  by  Yoffe  [41],  The  mechanism  of 
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Fig.  18  Crack  branching  after  a  mode  I  crack  initiation  if  the  notch  radius  was  small  but  impact 
loading  was  high 


crack  branching  subjected  to  a  high  loading  (impact)  rate  can  be  explained  using  the 
dynamic  energy  release  rate  (driving  force).  The  availability  of  kinetic  energy  due 
to  high  impact  speed  tends  to  create  more  fracture  surfaces  for  absorbing  energy. 
Therefore,  branched  cracks  easily  occur  in  a  high  impact  loading  case. 


3. 7  Dynamic  Crack  Kinking  and  Penetration  at  an  Interface 

3.7.1  Weak  Interfaces  with  Different  Interfaciai  Angies 

We  define  a  weak  interface  if  the  strength  or  toughness  of  this  interface  is  much  less 
than  that  of  the  strength  or  toughness  of  the  bulk  material  (here  is  Homalite).  The 
dynamic  crack  kinking  phenomenon  as  shown  in  Fig.  16  will  be  slightly  different  if 
the  interface  crack  is  changed.  For  the  case  of  the  interfaciai  angles  of  45°  [30],  the 
incident  crack  reached  the  interface  around  153  p.s  after  impact  as  seen  in  Fig.  19a. 

It  is  observed  that  the  symmetric  fringe  pattern  of  the  incident  mode-I  crack 
disappeared  as  soon  as  the  crack  deflected  into  the  interface.  The  shape  of  the 
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Fig.  19  Crack  deflection  process  at  a  weak  interface  (interface  angle  45°) 


fringe  pattern  of  Fig.  19b  suggests  that  this  interfacial  crack  is  shear-dominated  at 
the  latter  propagation  stage.  In  Fig.  19c,  it  is  interesting  to  observe  that  after  this 
shear-dominated  crack  propagated  some  distance  along  the  interface,  some  sec¬ 
ondary  cracks  formed  at  one  side  of  the  interface.  These  secondary  cracks  are  locally 
mode-I  and  form  on  the  tension  side  of  the  sheared  interface. 

They  form  after  the  dominant  crack  has  propagated  along  the  interface  and  thus 
after  the  interface  has  already  failed  in  shear.  These  types  of  secondary  cracks  that 
are  a  by-product  of  the  shear  crack  growth  along  interfaces  have  already  been  ob¬ 
served  experimentally  and  are  always  associated  with  dynamic  shear-dominated 
crack  growth  along  weak  interfaces.  As  the  interfacial  angle  is  changed  to  60°,  the 
dynamic  crack  deflection  behavior  is  slightly  altered,  as  seen  in  Fig.  20  on  the  crack 
speed  records. 

First,  we  notice  the  crack  speed  jump  across  the  interface  at  about  150  ps.  Obvi¬ 
ously,  the  initial  interfacial  crack  speed  of  700  m/s  is  much  higher  than  the  incident 
crack  speed  which  is  about  400  m/s  for  the  interfacial  angle  60°  case.  However,  the 
interfacial  crack  speed  reduced  to  350  m/s  soon  after  the  interfacial  crack  kinked 
into  the  right  side  of  monolithic  Homalite  (Fig.  16d).  The  experiment  also  suggests 
that  just  before  the  crack  kinking,  there  was  a  brief  crack  speed  reduction  charac¬ 
teristically  seen  in  several  failure  mode  transition  experiments  [45].  In  Fig.  20,  the 
comparison  of  the  crack  peed  history  is  made  for  two  different  interfacial  angle 
cases.  It  is  noticed  that  the  interfacial  crack  speed  for  the  interfacial  angle  45°  case 
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Fig.  20  Crack  speed  history  before  and  after  crack  deflection  at  a  weak  interface  (for  interfacial 
angles  45°  and  60°) 


is  always  higher  than  the  crack  speed  for  the  interfacial  angle  60°  case.  This  ex¬ 
perimental  phenomenon  was  analyzed  by  the  proposed  dynamic  fracture  mechanics 
model  [30], 


3.7.2  Modeling  of  Dynamic  Failure  Modes  Across  an  Interface 


For  the  incident  crack,  we  may  use  dynamic  fracture  mechanics  theory  [40]  to  fit 
the  stress  intensity  factors  A'/,  Kf,  and  the  T  stress  of  the  incident  crack.  Then,  we 
can  predict  the  key  parameters  of  the  kinked  interfacial  crack  using  fracture  me¬ 
chanics  theory  [47],  Our  first  step  is  to  obtain  the  static  counterparts  of  the  dynamic 
stress  intensity  factors.  After  crack  deflection,  the  static  stress  intensity  factors  of 
the  kinked  crack  can  be  calculated  using  these  relationships.  Let  Kf,  Ksnk  denote 
static  mode-I  and  mode-II  stress  intensity  factors  for  the  deflected  (kinked)  mixed¬ 
mode  crack,  and  they  are  related  to  the  static  stress  intensity  factors  of  the  incident 
dynamic  cracks  as  a  function  of  the  kinking  angle  f)  [1,48]: 


Kf  =  cuKs,+c]2Ki 


Ksk  =c2xK)+c22Ksu 


where  the  coefficients  are: 
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Let  K^k,  Kfk  be  dynamic  mode-I  and  mode-II  stress  intensity  factors  for  the 
deflected  (kinked)  mixed-mode  crack.  We  still  assume  that  the  universal  relation 
between  the  dynamic  and  static  stress  intensity  factors  also  holds  for  the  deflected 
crack  tip,  i.e., 

=  kuMKjf 

where  kj(v2)  and  kn(v2)  are  universal  functions  of  the  crack  tip  speed  defined  by 
Freund  [40].  Based  on  the  above  relations,  if  we  know  the  dynamic  stress  intensity 
factors,  the  crack  tip  speed  of  the  incident  crack,  the  kinking  angle  as  well  as  the 
crack  tip  speed  of  the  kinked  crack,  we  can  get  the  dynamic  stress  intensity  fac¬ 
tors  of  the  kinked  crack  and  hence  the  crack  tip  stress  fields  around  the  deflected 
crack  using  Eq.  7.  Then,  the  fringe  patterns  of  the  interfacial  crack  at  the  moment 
of  crack  deflection  can  be  predicted  using  Eqs.  1  and  5.  Because  it  is  very  hard  to 
record  the  exact  moment  of  crack  kinking  at  the  interface  in  dynamic  fracture  ex¬ 
periments,  average  values  of  the  stress  intensity  factors  of  the  incident  crack  were 
used  to  calculate  the  stress  intensity  factors  of  the  kinked  interfacial  crack. 

Figure  21  shows  the  predicted  fringe  pattern  of  a  kinked  interfacial  crack  (inter¬ 
facial  angle  60°).  The  coordinate  origin  is  located  at  the  intersection  point  of  the 
incident  crack  and  the  kinked  crack,  and  its  x-axis  is  along  the  interface.  For  this 
case,  the  crack  tip  speed  of  the  incident  crack  was  around  400  m/s. 

Right  after  crack  kinking,  the  interfacial  crack  tip  speed  was  about  700  m/s  [30]. 
Since  the  T-stress  of  the  incident  crack  is  around  1  MPa,  and  there  are  no  accurate 
results  for  the  T-stress  of  the  kinked  crack  [49],  the  T-stress  of  the  kinked  crack 
was  assumed  to  be  zero  in  all  our  predictions.  In  order  to  highlight  our  comparison 
of  the  predicted  and  experimental  fringe  patterns,  only  fringe  order  1  was  plotted. 
Figure  21a  presents  the  experimental  fringe  showing  the  transition  from  an  incident 
crack  to  an  interface  crack  (the  horizontal  line  was  the  interface).  Figure  21b  showed 
the  predicted  fringe  and  the  two  kinds  of  photoelasticity  fringes  were  very  similar  so 
our  dynamic  fracture  mechanics  modeling  and  assumptions  were  reasonable.  How¬ 
ever,  some  discrepancy  is  also  noticed  because  it  is  very  hard  to  take  one  photo  at 
the  right  time  and  right  position  as  the  theoretically  predicted  one. 

One  interesting  observation  is  the  large  concave  wedge  effect.  As  seen  in 
Fig.  21a,  two  caustic  spots  (one  caused  by  the  kinked  crack  tip  and  the  other 
due  to  a  large  concave  wedge)  were  clearly  observed  when  a  mode-I  incident  crack 
kinked  along  a  weak  interface  with  a  large  kinking  angle  (60°,  see  the  same  photo 
in  Fig.  16c). 

In  most  previous  crack  kinking  analyses,  researchers  only  considered  the  singu¬ 
lar  stress  field  due  to  a  kinked  daughter  crack  and  ignored  the  singular  stress  field 
of  a  concave  wedge.  Interestingly,  William's  classical  solution  of  wedge  stress  sin¬ 
gularities  [22]  is  the  foundation  of  the  full-field  stress  field  of  a  traction-free  crack 
in  Linear  Elastic  Fracture  Mechanics.  Indeed,  Cotterell  and  Rice’s  classical  work 
mainly  deals  with  a  slightly  kinked  crack,  not  a  large  kinking  angle  case  [50],  To 
authors’  knowledge,  only  Azhdari  and  Nemat-Nasser  provided  a  simple  explanation 
to  this  phenomenon  for  a  static  crack  kinking  case  [51]. 
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Fig.  21  Comparison  of  (a) 
experimental  fringe  (strong 
interface)  and  (b)  predicted 
fringe  pattern  (vj  =  400  m/s, 
v2  =  700  m/s,  N  =  1, 
interfacial  angle  60°) 


a 

Caustics  caused  by  the  Caustics  caused  by  the 

large  concave  wedge  interfacial  crack  tip 


3.7.3  Mode  Mixity  of  the  Kinked  Interfacial  Crack 

Mode  mixity  is  an  important  parameter  in  interface  fracture  analysis.  It  is  the  ra¬ 
tio  measure  of  the  mode  II  stress  intensity  factor/energy  release  rate  over  its  mode 
I  counterparts.  There  were  numerous  studies  on  mode  mixity  in  static  fracture  cases 
but  very  few  results  were  reported  in  dynamic  fracture  investigation  [52,53].  In¬ 
deed,  mode  mixity  is  a  key  parameter  in  controlling  failure  mode  transitions  along 
interfaces.  In  this  investigation,  when  the  incident  mode  I  crack  reached  the  inter¬ 
face,  it  kinked  along  the  interface  and  became  one  mixed  mode  crack.  Based  on 
Eq.  7,  for  the  kinked  interfacial  crack,  its  mode  mixity  depends  on  the  kinked  crack 
tip  speed  and  the  interfacial  or  kinking  angle.  Obviously,  the  dependence  on  the 
dynamic  mode  mixity  on  the  crack  tip  speed  is  a  special  phenomenon  in  dynamic 
fracture  mechanics.  The  variations  of  the  mode  mixity  with  the  interfacial  angle  and 
the  kinked  crack  tip  speed  are  plotted  in  Fig.  22.  It  is  not  surprising  to  see  that  when 
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Fig.  22  Dynamic  fracture  mode  mixity  as  functions  of  (a)  interfacial  angle  and  (b)  crack  tip  speed 


the  crack  tip  speed  of  the  interfacial  crack  remained  constant,  mode  mixity  increased 
with  the  increase  of  the  kinking  angle.  In  other  words,  the  larger  of  the  interfacial 
angle,  the  larger  is  the  mode  II  component  for  the  mixed  mode  crack.  This  result  is 
similar  to  the  common  conclusion  in  static  crack  kinking  analysis  [48]. 

As  a  special  feature  of  dynamic  crack  kinking,  the  mode  mixity  increases  w'ith 
the  increasing  kinked  crack  speed  if  the  interface  angle  is  fixed.  In  Fig.  22a,  for  a 
fixed  kinking  or  interfacial  angle  50°,  the  mode  mixity  for  a  high  crack  tip  speed 
(90%  of  the  shear  wave  speed  of  the  matrix  material)  is  almost  doubled  compared 
to  tlie  mode  mixity  for  a  static  kinked  crack.  Here,  we  should  notice  that  the  crack 
tip  speed  of  the  kinked  crack  is  related  to  the  interfacial  bonding  strength  [30].  A 
weak  interface  will  lead  to  a  fast  interfacial  crack  tip  and  a  high  mode-II  component 
as  a  result.  Figure  22b  shows  the  mode  mixity  dependence  on  the  crack  tip  speed 
for  different  kinking  angles.  It  is  interesting  to  see  that  each  curve  has  a  similar 
shape  and  is  shifted  by  some  amount  for  a  different  kinking  angle.  In  this  investi¬ 
gation,  the  kinking  or  interfacial  angle  is  limited  to  0-90°.  Recently,  Rousseau  and 
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Rosakis  [54]  examined  important  crack  kinking  behavior  for  very  large  interfacial 
angle  (greater  than  90°).  One  important  difference  is  that  their  incident  crack  was 
an  inter-sonic  shear  crack  along  a  weak  path  rather  than  a  slow  mode-I  crack  as  in 
our  investigation  In  order  to  suppress  possible  crack  branching,  our  incident  crack 
speed  was  controlled  to  be  less  than  the  crack  branching  speed  (around  30-40%  of 
the  shear  wave  speed  for  Homailte-100  [55]. 


3.8  Two-layer  Specimens  with  Direct  Impact  on  the  Brittle 
Polymeric  Layer  (Fig.  23) 

The  above  discussions  are  focused  on  the  three-layer  specimen  or  the  model  sand¬ 
wich  structure.  In  order  to  understand  the  dynamic  response  and  failure  of  general 
layered  materials  or  structures,  experimental  investigation  on  the  two-layer  speci¬ 
men  is  necessary  [56]. 

In  some  applications  of  layered  materials,  direct  impact  on  the  brittle  layer  such 
as  ceramics  layer  in  composite  armor  is  also  very  important  since  their  failure  be¬ 
havior  might  be  very  different  [57],  As  shown  in  Fig.  23  a  complicated  stress  wave 
propagation,  projectile  penetration  and  large  deformation  process  was  observed  for 
a  two-layer  system  with  a  weak  and  ductile  5083  bond.  As  soon  as  the  projectile 
impacted  the  transparent  Homalite  layer,  a  series  of  fringe  patterns  related  to  stress 
waves  developed  around  the  impact  site.  The  projectile  head  kept  moving  and  com¬ 
plicated  stress  wave  interaction  and  propagation  were  observed  around  10  [is  after 
impact  (Fig.  23b,  c). 

It  should  be  noticed  that  the  movement  of  the  projectile  was  not  obvious.  In  the 
next  stage  of  projectile  penetration,  a  black  contact  zone  was  observed  at  the  im¬ 
pact  site  and  an  intra-layer  crack  initiated  and  propagated  towards  the  interface.  In 
Fig.  23e,  it  seemed  that  this  intra-layer  crack  arrested  at  the  interface  since  no  signif¬ 
icant  interfacial  crack  propagation  (moving  fringe  concentration)  was  observed  in 
the  following  images.  As  discussed  by  Xu  and  Rosakis  [33],  the  crack  arrest  mech¬ 
anism  by  a  ductile  and  weak  adhesive  layer  is  due  to  the  dramatically  stress  wave 
gradient  change  across  the  adhesive  layer.  After  300  [is  of  impact  (Fig.  23f),  a  large 
contact  zone  appeared  in  a  growing  black  zone  connected  to  the  projectile.  Due 
to  large  deformation  inside  the  contact  zone,  those  transmitted  laser  rays  were  de¬ 
flected  and  cannot  enter  the  camera  so  only  a  black  zone  was  recorded.  Meanwhile, 
some  materials  were  flying  out  from  the  Homalite  layer  because  of  the  fragments 
of  brittle  Homalite  subjected  to  high  contact  force.  Compared  to  the  impact  on  the 
strong  steel  layer  of  previous  cases,  direct  impact  on  the  brittle  Homalite  layer  has 
more  severe  damage  in  the  local  impact  site. 
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Fig.  23  Failure  sequence  of  a  two-layer  specimen  (a)  impacted  directly  on  the  brittle  Homalite 
layer  (b-c)  show  a  complicated  stress  wave  process,  (d-f)  reveal  intra-layer  crack  initiation  and 
propagation  and  a  penetration  procedure  with  large  local  deformation 


4  Conclusions 

Inter-layer  crack  growth  (delamination)  is  the  dominant  dynamic  failure  mode  for 
layered  materials  and  sandwich  structures  subjected  to  out-of-plane  impact.  These 
cracks  appear  to  be  shear-dominated  and  proceed  with  intersonic  speeds.  Intra-layer 
cracking  always  occurs  soon  after  some  amount  of  inter-layer  delamination  has 
already  happened  and  proceeds  through  the  spreading  and  branching  of  local  mode 
I  cracks  into  the  core  layer.  Intra-layer  or  core/matrix  cracking  often  initiates  at 
the  interface  as  a  result  of  inter-layer  crack  kinking  into  the  core.  If  the  speed  of 
the  kinked  intra-layer  crack  reaches  a  critical  value,  multiple  crack  branching  may 
also  occur  inside  the  brittle  core  layer.  If  free-edge  effects  at  the  bimaterial  comers 
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are  eliminated  (long  model  sandwich  specimens),  the  failure  sequence  is  slightly 
modified.  Specifically,  the  inter-layer  cracks  initiate  from  positions  where  the 
inter-layer  shear  stress  reaches  a  local  maximum  equal  to  the  shear  strength  of  the 
bond.  These  cracks  create  a  local  shear-driven  delamination  directly  above  the  point 
of  impact.  Intra-layer  cracks  following  this  process  also  kink  from  these  positions 
into  the  sandwich  core. 

For  some  individual  dynamic  failure  modes,  we  find  that  the  dynamic  crack 
branching  can  exceed  the  limit  angles  45°  (theoretical  result)  for  a  second  dynamic 
crack,  and  the  crack  tip  speeds  before  and  after  the  crack  branching  are  almost  same. 
Crack  kink  and  penetration  at  interfaces  are  mainly  governed  by  the  ratio  competi¬ 
tion  of  the  dynamic  energy  release  rates  of  the  incident  crack  and  the  kinked  crack, 
and  the  dynamic  fracture  toughness  ratio  of  the  matrix  (bulk  material)  and  the  inter¬ 
face.  For  the  kinked  crack  speed,  it  increases  if  the  interfacial  angle  and  interfacial 
strength  decrease.  A  strength-based  criterion  is  proposed  to  predict  the  interfacial 
debonding  ahead  of  a  main  crack  w'ell. 
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Attending  in-class  paper  demonstrations  and  performing  experiments  in  solid  mechanics  courses 
ore  very  effective  ways  for  students  ta  gain  an  understanding  af  the  complicated  concepts  af 
mechanics.  This  paper  explores  a  few  applications  of  dogbane  tensile  tests,  stress  concentrations 
and  crack  kinking  ar  mixed-mode  fracture  Furthermore .  this  handy  technique  can  be  extended  to 
other  broader  areas  af  mechanics  education.  Since  only  simple  materials  and  supplies  are  used: 
copy  paper,  staples,  scissors  and  a  paper  punching  machine,  students  can  repeat  these  typical 
mechanics  experiments  in  future  in  other  locations ,  such  as  in  an  o  ffice  ar  at  home  Therefore .  this 
simple  and  effective  technique  can  have  a  remarkable  influence  an  the  student's  long-term  career 

Keywords:  mechanics  education;  mechanics  of  materials;  solid  mechanics,  experimental 
mechanics 


INTRODUCTION 

THE  AUTHOR’S  first-year  teaching  assignments 
at  Vanderbilt  University  included  two  graduate 
courses:  Advanced  Mechanics  of  Solids,  Parts  I 
and  II  Part  I  mainly  covers  the  theory  of  elasticity, 
w'hile  Part  II  ranges  from  plasticity,  buckling,  and 
stress  waves  to  fracture  mechanics.  These  highly 
theoretical  courses  require  a  strong  mathematical 
background  and  the  author  found  that  many 
graduate  students  felt  overwhelmed  by  the  deriva¬ 
tion  of  equations.  Therefore,  the  author  decided  to 
add  some  experimental  demonstrations  to  enable 
the  students  to  understand  complex  mechanics 
phenomena  Before  teaching  the  Rayleigh  surface 
wave,  the  author  showed  a  vivid  movie  on  stress 
wave  propagation  (which  can  now  be  seen  on  the 
author’s  Website)  based  on  the  author’s  previous 
unpact  mechanics  research  using  high-speed  photo¬ 
graphy  [I].  The  movie  captured  the  students'  inter¬ 
est,  hence  they  concentrated  on  the  long  derivation 
involving  the  speed  relationship  between  the 
Rayleigh  wave  and  the  shear  wave.  It  was  then 
that  the  author  recalled  a  famous  maxim,  “I  hear 
and  I  forget;  I  see  and  I  remember;  I  do  and  I 
understand,”  (Chinese  philosopher  and  educator 
Lao  Tsu,  604-531  BC).  Therefore,  he  decided  to 
develop  a  few  in-class  demonstrations,  both  for  the 
two  graduate  courses  mentioned  above,  and  for 
two  other  undergraduate  courses:  Mechanics  of 
Materials  and  Statics. 

CURRENT  CHALLENGE  OF  SOLID 
MECHANICS  COURSES 

Since  traditional  solid  mechanics  courses  involve 
the  successful  applications  of  applied  mathematics 
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[2-3],  they  arc  challenging  for  the  engineering 
students  who  have  less  mathematical  training. 
For  example,  students  who  have  never  taken  a 
Partial  Differential  Equation  course  would  find 
that  it  is  difficult  to  understand  Boundary  Value 
Problems  (BVP,  a  pure  mathematical  terminology') 
in  the  theory  of  elasticity  [4-5]  Therefore,  for 
modern  mechanics  educators,  the  use  of  too 
much  mathematical  derivation  without  the  provi¬ 
sion  of  engaging  examples  may  turn  the  students' 
interests  to  other  fields  that  they  perceive  to  be 
more  fun;  this  docs  not  help  the  development  of 
mechanics  disciplines.  Fortunately,  new  informa¬ 
tion  technology  provides  us  with  excellent  tools  for 
dynamic  instruction  and  effective  visualization. 

The  main  objective  of  this  paper  is  to  demon¬ 
strate  that  a  solid  mechanics  course  is  a  series  of 
dynamic  and  fun  events,  rather  than  a  cluster  of 
bonng  equations.  Indeed,  in  his  series  of  classical 
mechanics  textbooks,  outstanding  mechanics 
researcher  and  educator,  Professor  S.  Timoshenko, 
often  used  experimental  stress  techniques  such  as 
photoelasticity  to  explain  complicated  solid 
mechanics  problems  [4].  Now,  modern  multimedia 
could  be  used  to  realize  this  goal  more  effectively. 
Usually,  we  want  to  expose  students  to  three 
areas  of  mechanics — experimentation,  theory  and 
computation  [6].  The  use  of  multimedia,  such  as 
personal  computers  and  the  Internet,  can  be  very 
helpful  in  explaining  theoretical  and  computa¬ 
tional  aspects  of  mechanics,  and  can  even  be 
used  to  conduct  “virtual  experiments"  [7]. 
However,  for  undergraduate  introductory  level 
classes  such  as  the  mechanics  of  materials,  the 
impact  of  multimedia  on  teaching  fundamental 
concepts  has  been  less  successful  [8].  Therefore, 
effective  and  original  experimental  education  such 
as  doing  tensile  tests  in  mechanics  of  material 
courses  still  cannot  be  replaced  [9-1 1]. 
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IN-CLASS  PAPER  DEMONSTRATIONS 
AND  EXPERIMENTS 

For  experimental  demonstrations,  not  every 
educator  has  sufficient  laboratory  facilities  or 
experience  to  conduct  educational  experiments 
such  as  the  classical  photoelasticity  technique  [9]. 
In  order  to  develop  a  handy  education  technique, 
the  author  has  explored  an  in-class  paper  (A4  copy 
paper,  etc.)  experiment,  which  can  be  employed 
not  only  in  classroom  teaching,  but  also  for 
outreach  presentations.  Since  the  mechanics  of 
paper  is  a  typical  plane  stress  problem,  the 
author  has  extensively  used  A4  copy  paper  in 
classroom  lectures  for  mechanics  of  materials 
(dogbone  tensile  tests),  theory  of  elasticity  (stress 
concentration  around  a  central  hole),  plates  and 
shells  (KirchhofTs  assumptions)  and  fracture 
mechanics  (crack  kinking  and  fracture  modes). 
For  example,  a  classical  mixed-mode  fracture 
experiment  [12]  can  easily  be  conducted  and 
takes  only  a  few  minutes.  For  a  laboratory  experi¬ 
ment  to  demonstrate  the  same  mixed-mode  frac¬ 
ture  phenomenon,  the  process  may  take  at  least 
several  days,  from  specimen  preparation  to  the 
final  experiment.  Moreover,  the  cost  would  be 
much  higher  than  that  of  our  paper  experiment. 
In  the  following  sections,  we  will  describe  in  detail 
the  steps  of  a  few  typical  in-class  experiments  from 
specimen  preparations  to  final  experiments,  and 
give  the  related  theoretical  background. 

Applications  in  basic  mechanics  courses — dogbone 
tensile  tests  and  stress  concentrations 

Traditional  mechanics  of  materials  courses 
mainly  cover  the  mechanical  response  of  a  metal 
bar  subjected  to  tension  and  other  loadings.  There¬ 
fore,  the  corresponding  laboratory  class  is  focused 
on  the  stress  and  strain  measurements  of  a  steel  or 
aluminum  cylinder  specimen.  With  the  extensive 


applications  of  new  polymers  and  composite 
materials,  as  well  as  thin  films  [13],  dogbone 
specimens  have  become  quite  popular  for  use  in 
tensile  experiments.  Figure  1  shows  the  prepara¬ 
tion  process  of  a  paper  dogbone  specimen.  The 
first  step  is  to  print  out  symmetrical  axes  and 
specimen  edges  using  computer  software  such  as 
Microsoft  Word.  The  second  step  is  to  fold  the 
whole  paper  into  a  quarter  so  that  only  “Part  1” 
can  be  seen.  Then,  using  scissors,  the  curved  edges 
are  carefully  cut.  Now  a  dogbone  specimen  is 
ready  to  be  tested  as  seen  in  Fig.  1(c).  For  compar¬ 
ison,  a  rectangular  paper  specimen  (a  strip)  is  also 
cut.  Students  can  pull  the  ends  of  these  two  kinds 
of  tensile  specimens  and  compare  their  failure 
patterns.  The  final  fracture  of  the  dogbone  speci¬ 
men  is  often  at  the  specimen  center  (a  valid  tensile 
test).  However,  the  location  of  the  final  fracture  of 
the  rectangular  paper  specimen  is  random,  i.e.,  at 
the  specimen  center  or  the  grip  location,  etc. 

The  dogbone  specimen  can  be  used  to  demon¬ 
strate  another  important  solid  mechanics  phenom¬ 
enon:  stress  concentration  As  shown  in  Fig.  2(a),  a 
paper  punch  is  used  to  punch  a  hole  in  the 
dogbone  specimen.  According  to  classical  elasticity 
solutions  for  an  infinite  plate  with  a  circular  hole 
of  radius,  a,  subjected  to  remote  uniaxial  tension, 
cr,  the  stress  components  along  the  hole  edge 
expressed  in  a  polar  coordinate  system  are  [4] 


*"  =  ?( 1  -3  )  +?( i  +7r-7r)cos(20) 

Vm  =^t  (  1  +3  )  --z  (  1  +~x)  cos(20)  (1) 


For  6  =  7r/2,  37t/2,  ct^  attains  its  maximum  value 
3cr.  Here,  “3”  is  called  the  stress  concentration 


(a)  (b)  (c) 

Fig.  1 .  Dogbone  tension  specimen  preparation  using  A4  copy  paper:  (a)  mark  fold  lines  on  A4  paper  using  computer  software  such  as 
Microsoft  Word;  (b)  fold  the  whole  paper  into  one  quarter  and  cut  paper  using  scissors:  (c)  a  complete  dogbone  specimen. 
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factor.  For  a  finite  dogbone  specimen  with  a 
central  hole,  the  stress  concentration  factor  is  less 
than  3.  By  conducting  tensile  tests  of  the  dogbone 
specimens  with  and  without  central  holes,  students 
will  be  able  to  notice  that  the  dogbone  specimen 
with  a  central  hole  often  fails  at  the  hole  edge  due 
to  the  high  stress  concentration.  Because  these 
elasticity  equations  are  usually  beyond  the  under¬ 
graduate  student  curriculum,  the  use  of  the  stress 
flow  concept  is  more  effective  to  explain  the  stress 
concentration  phenomenon.  As  with  fluid  visual¬ 
ization,  the  tensile  stress  of  the  dogbone  specimen 
at  the  far  field  is  quite  uniform  as  shown  in  Fig 
2(b).  However,  the  stress  flow  becomes  “very 
crowded"  when  it  encounters  the  central  hole, 
since  the  stress  flow  cannot  pass  across  the  hole. 
Therefore,  a  stress  concentration  occurs  at  the 
location  with  a  significant  flow  concentration  as 
seen  in  Figure  2(b).  This  is  a  clear  explanation  of 
the  stress  concentration  for  undergraduate 
students  (even  for  graduate  students).  The  stress 
flow  concept  could  be  more  effective  if  students 
have  taken  the  fluid  mechanics  course,  and  know- 
how  to  apply  the  continuity  or  mass  equation 

Applications  in  advanced  mechanics  courses — 
fracture  and  fiber-reinforced  composites 

For  advanced  solid  mechanics  courses,  fracture 
mechanics  is  one  of  the  most  important  subjects. 
After  the  theory  of  elasticity  and  plasticity,  the 
major  challenge  for  graduate  students  is  to  under¬ 
stand  the  concepts  of  the  singular  stress  field  at  a 
crack  tip  (theoretical  stress  values  are  infinity, 
based  on  the  elasticity  solutions),  and  the  stress 
intensity  factor  of  a  crack  in  Linear  Elastic  Frac¬ 
ture  Mechanics.  In  experimental  mechanics,  the 
singular  stress  field  will  result  in  a  high  stress 
concentration  at  the  crack  tip  or  bi-material 
comer  [14].  Of  course,  a  fracture  experiment 
using  photoelasticity  or  other  techniques  may 


lead  to  a  much  clearer  understanding  of  fracture 
mechanics.  However,  many  instructors  do  not 
have  sufficient  time  and  laboratory  facilities  to 
prepare  these  demonstrations.  The  author  has 
developed  a  convenient  paper  experiment  to 
demonstrate  the  concepts  of  the  stress  intensity 
factor  and  mixed-mode  fracture,  which  are  among 
the  most  important  topics  in  fracture  mechanics 
courses.  As  shown  in  Fig.  3,  first  we  fold  a  piece  of 
A4  copy  paper  and  cut  an  initial  ‘crack',  then  we 
load  the  initial  ‘crack'  by  pulling  on  the  paper. 
When  the  pulling  load  is  sufficiently  large,  sudden 
crack  kinking  occurs,  due  to  the  mixed-mode  load, 
as  shown  in  Frig.  3(c).  This  experiment  is  easy  to 
conduct  and  one  can  change  the  fold  angle  of  the 
initial  crack,  and  obtain  different  crack  kinking 
patterns.  Indeed,  the  crack  kinking  phenomenon 
involves  an  interesting  theoretical  background. 

A  crack  will  propagate  away  from  its  original 
path  due  to  any  change  in  local  fracture  mode 
mixity  (a  mechanics  factor)  or  fracture  toughness 
(a  material  factor).  This  phenomenon  is  often 
called  “crack  kinking  or  deflection"  [15-18].  For 
a  two-dimensional  elastic  solid,  the  full-field  stress 
tensor  of  a  mixed-mode  main  crack  as  seen  in  Frig. 
3(b),  can  be  expressed  in  a  polar  coordinate  system 
as  in  Williams  [19]  and  Anderson  [20]. 


(2) 


K  '  / 

ism  .. 

+  -7r=£i A9)  +  o(A)  (/,;=  1,2 

V  2nr  J 


w  here  m  and  k  denotes  “main  crack"  and  “kinked 
crack"  respectively.  K\  and  K\\  are  mode  I  and 
mode  II  stress  intensity  factors;  T  is  a  nonsingular 
term;  0(/4)  represents  higher  order  terms  of  the 
length  scale  r  and  will  be  dropped  if  the  kinked 
crack  length  /  is  very  small;  and  known  functions 
Ef  (0),  (0)  represent  the  angular  variations  of 


Strew 

concentrations 


Fig.  2.  Stress  concentrations  during  tension  experiments,  (a)  punch  a  hole  in  a  dogbone  specimen  using  a  paper  punching  machine:  (b) 

stress  How  around  the  circular  hole. 
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Fig.  3.  Paper  experiment  lo  demonstrate  lhc  crack  kinking  phenomenon;  (a)  folding  paper  and  cut  an  initial  ‘crack’,  (b)  loading  the 
initial  ‘crack'  by  pulling  paper;  (c)  crack  kinking  due  to  mixed-mode  loading. 


stress  components,  Previous  research  on  crack 
kinking  reveals  that  there  is  a  relation  between 
the  stress  intensity  factors  before  and  after  crack 
kinking,  which  can  be  expressed  in  a  combined 
way  as  [16,  21]: 

K^culCr  +  cnKJl  +  biTy/l  (3) 

til  =  cuKf  +  cl2K”  +  bi  Tsfl.  (4) 

Numerous  results  have  been  reported  to  determine 
the  coefficients  Cy  and  bt  ( ij  =  1,2).  From  the  above 
relations  it  is  obvious  that  the  stress  intensity 
factors  of  the  kinked  crack  depend  on  the  stress 
intensity  factors  of  the  initial  main  crack  Further¬ 
more,  we  can  use  stress  or  energy  criteria  to  predict 
the  crack  kinking  initiation  load  and  its  angle  of 
kink  [12,  22].  Therefore,  by  using  different  kinds  of 
paper  materials  and  differently  inclined  angles  of 


the  main  cracks,  a  course  project  on  mixed-mode 
fracture  allows  students  to  obtain  a  deep  under¬ 
standing  of  the  concepts  of  fracture  mechanics.  In 
addition  to  demonstrating  the  mixed-mode  frac¬ 
ture  phenomenon,  paper  with  an  edge  crack  as 
shown  in  Fig.  4(a)  can  be  used  to  demonstrate 
mode  I-opening,  mode-II  in-plane  shear  and 
mode-III  out-of-plane  shear,  by  holding  the  two 
edges  and  applying  different  kinks  of  load  (in¬ 
plane  tension;  tearing,  etc). 

The  edge-cracked  paper  could  be  slightly  modi¬ 
fied  to  demonstrate  the  mechanics  of  fiber-rein- 
forced  composites.  As  seen  in  Fig.  4(a),  several 
staples  fixed  in  front  of  a  main  crack  are  used  to 
increase  the  local  fracture  resistance  of  the  paper 
material.  These  stiff  and  strong  staples  are  analo¬ 
gous  to  short  micro/nano  fibers  used  in  composites 
[23-24].  In  Fig.  4(b),  after  we  pull  the  paper  with  an 
edge  crack,  the  main  crack  kinks  along  the  ends  of 
the  short  staples,  since  the  local  fracture  toughness 


MODE-1 


t 


(a) 


(b) 


Fig.  4.  Demonstration  of  fracture  modes  and  fiber-reinforced  composites  using  paper  and  siaples.  (a)  mode-I  opening  and  mode-II 
in-planc  shear  modes;  (b)  short  fibers  can  provide  high  fraciurc  resisiance. 
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of  the  paper  is  significantly  increased  by  these  strong 
staples.  This  experiment  demonstrates  a  fracture 
mechanics  conclusion;  a  crack  tends  to  choose  a 
path  of  low  fracture  resistance  (toughness). 

These  paper  experiments  and  demonstrations 
used  in  the  classroom  may  have  a  long-term 
influence  on  the  student's  future  career.  Since  we 
need  only  simple  materials  and  supplies;  copy 
paper,  staples,  scissors  and  paper  punching 
machines,  students  can  repeat  these  experiments 
elsewhere  in  the  future,  such  as  in  an  office  or  their 
homes.  Therefore,  the  general  public  can  easily 
understand  the  important  role  of  solid  mechanics 
in  scientific  research. 


CONCLUSIONS 

ln-class  paper  experiments  for  solid  mechanics 
courses  are  an  effective  way  for  students  to  learn  to 
understand  complicated  mechanics  concepts  This 
paper  explores  a  few  such  applications  in  dogbone 
tensile  tests,  stress  concentrations  and  crack  kink¬ 
ing  or  mixed-mode  fracture.  This  handy  technique 
of  solid  mechanics  education  can  be  used  to 
demonstrate  the  impact  of  solid  mechanics 
research  to  the  general  public. 
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Abstract  An  integrated  experimental  and  numerical  anal¬ 
ysis  is  carried  out  to  study  the  interfacial  shear  strength  of 
bonded  materials.  Two  types  of  shear  tests,  namely  the 
losipescu  shear  test,  and  the  short-beam  shear  test  are 
employed  to  understand  the  effect  of  interfacial  stress  on 
the  interfacial  shear  strength  measurements.  The  measured 
average  shear  strengths  are  very  close,  even  though  the 
interfacial  shear  stress  distributions  of  these  two  kinds  of 
specimens  are  very  different.  Therefore,  we  conclude  that 
the  interfacial  stress  distribution  has  the  least  effect  on  the 
interfacial  strength  measurement  if  the  interfacial  shear 
stress  is  non-singular. 

Keywords  Adhesion  •  Interfacial  strength  Photo-elasticity 
Finite  element  simulations  *  Monte  Carlo  simulations  * 

Shear  failure 

Introduction 

Composite  materials  find  extensive  applications  in  many 
engineering  fields  [1].  Since  composites  are  made  of  two  or 
more  materials,  there  exists  an  interface  between  the 
different  types  of  materials.  This  interface  is  usually  the 
weakest  link  and  hence  it  becomes  very  important  to 
characterize  its  mechanical  properties  [2].  For  example,  the 
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interfacial  strength  is  an  important  parameter  which 
determines  the  interface  failure  subject  to  dynamic  loading 
[3—5]  Recent  atomistic  simulations  show  that  the  interfacial 
shear  strength  of  the  copper-niobium  interface  is  always 
lower  than  the  shear  strength  of  perfect  crystals,  and  is 
spatially  non-uniform  [6]  There  are  several  different 
approaches  to  measure  the  interfacial  shear  strength  of 
composites  and  bonded  materials  such  as  the  fiber  push-out 
test  for  the  fiber  'matrix  interface  [7,  8].  Recently,  the  issues 
on  the  free-edge  stress  singularity  and  the  highly  non- 
uniform  stress  distribution  have  been  raised  for  this  type  of 
an  experiment  [9].  An  important  question  arising  out  of  this 
discussion  would  be  to  determine  if  the  interfacial  stress 
distribution  has  an  influence  on  the  interfacial  strength 
measurement.  In  order  to  get  nd  of  the  bimaterial  free-edge 
stress  singularity  and  address  this  fundamental  issue,  here  a 
same  material  joint  with  a  weakening  interface  is  used 
intentionally  This  ty  pe  of  material  joint  was  proposed  to 
investigate  the  dynamic  fracture  [10,  11].  In  case  of  shear 
specimens,  the  losipescu  shear  and  short-beam  shear 
specimens  will  be  employed  as  shown  in  Fig.  1 . 

The  losipescu  shear  test  provides  a  region  of  uniform 
shear  stress  in  the  test  zone  compared  to  a  parabolic  shear 
stress  variation  shown  by  the  short-beam  shear  specimen. 
The  losipescu  shear  test  was  developed  to  measure  the 
strength  of  composites  and  has  been  extensively  modified 
before  coming  to  its  present  state  [12-15].  Concurrently, 
the  short-beam  shear  test  is  also  used  widely  to  measure  the 
shear  strength  of  composite  materials.  The  losipescu 
specimen  shows  a  near  constant  variation  of  interfacial 
shear  stress,  which  might  be  seen  as  a  major  advantage  over 
other  types  of  shear  tests.  Therefore,  if  we  use  these  two 
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Fig.l  Variation  of  interfacial  shear  stress  for  (a)  Short-beam  shear 
specimen  (Parabolic  shear  stress  variation)  and  (b)  losipescu  shear 
specimen  (Near  constant  shear  stress  variation).  A  schematic  of  the 
boundary  conditions  used  in  finite  elements  (anti-symmetric)  is  also 
shown 


different  approaches  to  measure  the  interlaminar  shear 
strength  of  the  same  graphite/epoxy  material  system,  should 
we  get  very  different  shear  strengths  using  these  two 
approaches  due  to  very  different  interfaeial  shear  stress 
distributions,  or  almost  the  same  strength  because  the 
interfaeial  shear  stress  distribution  has  the  least  effect  on 
the  final  shear  strength? 

Surprisingly  no  one  has  explored  this  important  and 
fundamental  issue  over  the  past  40  years  (Adams  DF, 
Private  Communications,  2003).  Therefore,  in-depth  me¬ 
chanics  and  materials  work  are  much  needed  sinee  there  are 
several  other  approaches  to  measure  the  shear  strengths  of 
composite  materials  [1,  16].  The  objective  of  this  study  is  to 
analyze  the  effect  of  the  interfaeial  stress  distribution  on  the 
strength  measurements.  In  order  to  avoid  complicated 
failure  modes  during  shear  tests  [17],  and  to  only  focus 
on  the  shear  failure  of  an  interface,  we  design  and  test 
bonded  polymers  and  metals  such  that  the  bonded  interface 
is  the  only  weak  path  leading  to  the  final  shear  failure. 


Experimental  Investigation 

Test  specimens  were  made  of  three  different  types  of 
materials  including  Aluminum,  Polymethyl  methacrylate 


(PMMA)  and  Polycarbonate.  The  specimens  were  not 
monolithic  and  were  bonded  from  individual  halves. 
The  adhesive  used  in  this  experiment  is  Loctite  384. 
One  of  the  key  requirements  in  choosing  the  adhesive 
was  that  its  bonding  strength  should  be  far  lesser  than 
that  of  the  bulk  material.  The  specimens  were  bonded  at 
room  temperature  using  a  special  fixture  to  guarantee 
their  dimensionality.  The  individual  halves  were  sand¬ 
blasted  at  the  bonding  surfaces  in  order  to  provide  good 
adhesion. 

The  test  set-up  consists  of  three  parts  including  a 
mechanical  system  to  load  the  specimen,  an  optical  system 
to  develop  fringe  patterns  and  an  imaging  system  to  record 
the  images  [18].  The  mechanical  testing  system  included  an 
MTS  810  test  machine  and  an  losipescu  test  fixture.  The 
fringe  patterns  developed  only  in  ease  of  the  transparent 
Polycarbonate  specimens,  and  the  isoehromatie  fringe 
patterns  are  the  contours  of  the  maximum  in-plane  shear 
stress  [19]: 

_  (<7i  -  a2)  _  Nfa 
W  “  2  “  2  h 

where  crt  and  a2  arc  the  in-plane  principal  stresses,  N  is  the 
fringe  order,  fa  is  the  stress-fringe  constant,  and  h  is  the 
thickness  of  the  specimen  The  imaging  system  included  a 
high-resolution  digital  camera  to  capture  the  fringe  devel¬ 
opment  and  a  density’  filter  in  front  of  the  camera  to  reduce 
the  laser  intensity. 

Load  was  applied  to  the  movable  half  of  the  losipescu 
test  fixture  in  the  form  of  displacement  at  a  rate  of 
1  mm/min  until  failure.  About  25-35  specimens  of  each 
type  (losipescu  and  short-beam  shear)  were  tested  in 
order  to  ensure  repeatability.  Almost  all  of  the  specimens 
failed  due  to  shear  at  the  interface  which  was  necessary 
in  this  case  as  we  measure  the  shear  strengths.  Pictures 
of  the  fringe  patterns  at  load  levels  of  12.5  %  (TOO  N) 
and  25  %  (200  N)  of  failure  load  were  taken  for 
Polycarbonate  as  seen  in  Figs.  2  and  3,  and  videos 
depicting  the  development  of  fringes  from  the  beginning 
of  loading  till  failure  w'ere  recorded. 


Finite  Element  Modeling 


A  finite  element  model  of  the  specimen  wras  built  using  the 
software  ANSYS  11.0  to  obtain  the  stress  distribution  of 
the  bonded  specimen.  A  two  dimensional  analysis  was 
considered  for  bonded  Polycarbonate  materials  (Elastic 
Modulus  of  2.4  GPa,  and  a  Poisson's  Ratio  of  0.35).  A 
linear  elastic  and  isotropic  material  model  was  used  in  this 
analysis.  The  dimensions  of  the  specimen  were:  width  of 
76.2  mm  (individual  parts  have  a  width  of  38.1  mm)  and  a 
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Fig.  2  Comparison  between 
experimental  and  finite  element 
generated  fringe  patterns  for 
losipescu  specimen  made  of 
Polycarbonate 


Experimental 


Simulated 


(c)  Applied  Load  =  200  ±  5  N  (d) 


height  of  19.1  mm.  Plane  42  elements  were  exclusively 
used  in  order  to  meet  the  requirements  of  the  plotting 
software  Tccplot.  The  procedure  for  simulating  the 
boundary  conditions  was  adopted  similar  to  the  method 
reported  in  [17]. 

In  order  to  incorporate  a  realistic  simulation  of  the 
load  transfer  from  the  fixture  to  the  specimen,  an 
iterative  procedure  was  adopted.  Displacement  con¬ 
straints  are  applied,  similar  to  what  is  observed  in 
reality.  The  left  portion  of  the  specimen  is  fixed  and 
cannot  move  in  the  vertical  dimension.  The  top  left 
edge  node  is  restrained  from  moving  in  both  x  and  y 
directions  in  order  to  prevent  rigid  body  movements.  On 


the  right  portion,  an  initial  vertical  uniform  displace¬ 
ment  is  applied  Now  the  reaction  loads  at  each  of  the 
nodes  with  a  constraint  is  checked  to  verify  that  they 
are  not  in  tension  (in  real  experiments,  only  compres¬ 
sive  loading  is  obtained).  The  constraint  was  removed 
from  those  nodes  which  showed  a  tensile  reaction  force 
instead  of  a  compressive  one,  and  a  new  analysis  with 
the  updated  set  of  boundary'  conditions  was  carried  out. 
This  was  done  until  convergence  and  the  total  applied 
load  was  obtained  by  integrating  the  shear  stress  across 
the  bonded  interface  of  the  specimen.  This  procedure 
was  iteratively  carried  until  the  required  load  (for 
matching  the  fringe  patterns)  was  obtained.  It  should 
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Fig.  3  Comparison  between 
experimental  and  finite  element 
generated  fringe  patterns  for 
short-beam  shear  specimen 
made  of  Polycarbonate 
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be  noted  that  the  loading  is  anti-symmetric  and  remains 
so  after  every  iteration. 

The  numerical  photoelasticity  fringe  patterns  were 
obtained  by  plotting  the  stress  distribution  from  the 
finite  element  analysis  using  the  plotting  software 
Tecplot.  The  maximum  shear  stress  was  obtained  from 
the  finite  element  analysis  by  using  the  principal  stress 
values.  Using  equation  (1),  this  stress  value  was  then 
converted  to  a  fringe  order  N  which  in  turn  was  converted 
into  a  grayscale  value.  Half  order  fringes  (0.5,  1.5,  2.5, 
etc.)  were  given  a  value  of  255  and  full  order  fringes  (0,  1, 
2,  etc.)  were  assigned  a  value  of  0  on  the  grayscale 
spectrum. 


Results  and  Discussion 

A  direct  comparison  between  the  experimental  fringe 
patterns  and  the  fringe  patterns  is  obtained  from  the 
numerical  analysis  as  shown  by  the  results  presented  in 
Fig.  2  for  the  Iosipcscu  specimen  and  in  Fig.  3  for  the 
short-beam  shear  specimen.  At  lower  loading  levels 
(around  100  N),  the  fringes  begin  to  appear  from  the 
points  of  stress  concentration  (loading  point)  and  begin 
to  spread  out  towards  the  ends  of  the  shear  specimens. 
More  fringes  continue  to  appear  as  the  loading  is 
increased  and  the  patterns  arc  recorded  on  a  video 
camera  until  failure.  The  finite  element  analysis  captures 
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all  the  major  features  seen  in  the  experimental  patterns 
thereby  validating  our  finite  element  model  and  loading 
simulation. 

The  interfacial  shear  strengths  of  the  three  kinds  of 
bonded  materials  are  presented  in  Table  1.  It  was 
observed  that  the  average  shear  strengths  of  the  two 
types  of  different  shear  tests  w^ere  close.  In  the  case  of 
bonded  aluminum  specimens  it  is  seen  that  the  differ¬ 
ence  in  means  (between  Iosipescu  and  short-beam  shear 
specimens)  is  much  lesser  compared  to  the  other  two 
bonded  specimens.  In  general,  acrylic  materials  (PMMA) 
and  metals  show  a  better  bonding  in  comparison  with 
Polycarbonate  which  consistently  gives  a  lower  value  of 
the  interfacial  shear  strength.  Since  these  two  types  of 
shear  tests  lead  to  very  similar  results,  the  short-beam 
shear  test  will  have  advantages  like  low  machining  and 
fixture  cost,  and  shorter  time  for  preparation.  In  order  to 
provide  mechanics  insight  into  the  shear  failure  process 
at  the  interface,  Monte  Carlo  simulations  using  the 
initial  interfacial  flaw  assumption  were  carried  out 
[20].  The  interfacial  shear  stress  distributions  at  failure 
loads  of  the  Iosipescu  (with  same  cross-sectional  area  as 
the  short-beam  shear  specimen)  and  the  short-beam  shear 
specimens  from  the  finite  element  analysis  are  shown  in 
Fig.  4.  It  can  be  seen  that  the  variation  in  shear  stress  at 
the  interface  is  parabolic  for  the  short-beam  shear 
specimen,  but  almost  constant  for  the  Iosipescu  specimen. 
Therefore,  we  should  briefly  review  the  strength  concep¬ 
tion,  and  the  stress  distribution  of  the  specimens  to 
measure  the  strength. 

The  tensile  strength  definition  of  a  homogeneous 
material  is  quite  straightforward:  it  is  defined  as  the 
failure  stress  of  a  specimen  with  a  uniform  stress 
distribution  across  the  cross  section.  Since  the  final 
failure  of  the  material  is  associated  with  the  initial 
defect  distribution  (e.g.,  Weibull  strength  distribution)  at 
the  cross  section,  specimens  with  different  sizes  will 
yield  different  strength  data.  This  phenomenon  may  be 
defined  as  ‘the  material  size  effect.’  In  this  investiga¬ 
tion,  we  employed  homogeneous  materials  but  with 
artificial  interfaces  to  ensure  that  failure  will  occur 
along  these  weak  interfaces,  thereby  avoiding  compli¬ 
cated  bi-matcrial  problems  in  obtaining  useful  funda¬ 
mental  results  [21]. 


Fig.  4  Comparison  of  the  interface  shear  stress  distnbulions  at  failure 
load  of  728  N  between  the  short-beam  shear  specimen  and  the 
Iosipescu  shear  specimen  from  ihe  finile  element  analysis 

However,  for  bi-materials  specimens,  the  interfacial 
strength  measurement  becomes  quite  complicated  due  to 
existence  of  possible  free-edge  stress  singularity  [18]. 
Therefore,  it  is  necessary  to  examine  any  interfacial 
strength  experiment  carefully,  as  otherwise  it  might  yield 
meaningless  data  [9].  For  example,  the  mechanics  size 
effect  (related  to  the  stress  singularity)  might  become 
coupled  with  the  material  size  effect,  and  thus  lead  to 
complexities  in  interfacial  strength  measurement  data. 
Hence,  the  measured  nominal  interfacial  strengths 
cannot  be  used  in  mechanics  predictions  because  the 
interfacial  properties  obtained  from  laboratory  tests  arc 
very  different  from  the  real  values  of  structures  in 
service  (size-dependent).  If  the  free-edge  stress  singu¬ 
larity  were  to  be  successfully  removed  [21],  a  question 
would  be  whether  a  uniform  interface  stress  distribution 
is  required  to  measure  the  interfacial  strength.  The 
current  investigation  pointed  out  that  the  interfacial 
stress  distribution  has  the  least  effect  on  the  interfacial 
strength  measurement  as  long  as  the  stress  is  not 
singular  More  experimental  investigations  are  needed 
to  verify  whether  this  is  a  correct  statement  for  all  kinds 
of  materials  or  not.  The  current  work  provides  design 
guidance  for  future  development  of  all  interfacial 
strength  approaches. 


Table  1  Measured  interfacial 
bonding  strengths  (MPa) 


Specimen  Type 

Iosipescu  Shear 

Short-Beam  Shear 

Difference 

Aluminum-Aluminum 

10.75  ±2.39 

10  16±2.41 

5.5  % 

Polycarbonate-PC 

10.99±  1 .45 

8  51  ±1.13 

22.5  % 

PMMA-PMMA 

11.58±2.15 

10  19=0.57 

12% 
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Abstract 

In  this  paper,  an  efficient  method  to  characterize  the  mode  II  fracture  toughness,  Knc  of  materials 
with  preferred  interfaces  such  as  wood,  fiber  composite  materials  and  bonded  materials  based  on 
the  Iosipescu  fixture  is  developed.  Our  method  is  easy  to  implement  and  is  free  of  friction 
between  the  cracked  faces  thereby  offering  a  more  reasonable  value  of  Knc  than  previous 
methods.  Another  feature  of  our  new  approach  is  the  measurement  of  accurate  crack  initiation 
load.  Additionally,  a  calibration  chart  for  composite  materials  is  generated.  Finite  element 
analysis  along  with  cohesive  elements  is  used  to  validate  experimental  load-displacement  curves 
which  are  obtained  by  testing  bonded  polymers. 

1.  Introduction 

The  mode  II  fracture  toughness  has  long  been  a  difficult  property  to  quantify  for  isotropic 
materials  as  a  mode  II  crack  tends  to  kink  away  from  the  main  crack  (Anderson.  2004),  and 
thereby  makes  the  measured  fracture  toughness  indeed  the  mode  I  fracture  toughness.  A  mode  II 
crack  is  possible  only  for  special  materials  with  preferred  interfaces  like  wood,  bonded  materials 
or  composite  materials.  For  unidirectional  composite  materials  with  fiber/matrix  interfaces,  a 
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mode  II  crack  is  possible  due  to  the  strong  constraining  of  fibers  to  suppress  any  kinked  cracks, 
as  has  been  shown  in  recent  dynamic  shear  crack  experiments  by  Coker  and  Rosakis  (2001). 
Therefore,  several  experimental  approaches  have  been  proposed  to  measure  the  mode  II 
interlaminar  fracture  toughness  of  fiber-reinforced  composite  materials  such  as  the  end-notched 
flexure  test  (Daniel  and  Ishai.  2005).  But,  these  beam  bending  tests  usually  require  numerical 
procedures  like  data  reduction  methods  to  correct  for  the  shear  and  crack  tip  deflection 
(Yoshihara  and  Satoh.  2009).  These  issues  tend  to  complicate  the  experimental  measurement  of 
the  mode  II  fracture  toughness  which  is  a  fundamental  material  property. 

A  critical  issue  in  beam  bending  experiments  in  conjunction  with  Linear  Elastics  Fracture 
Mechanics  (LEFM)  analysis  is  the  significant  absence  of  the  “K-dominance  zoncv  in  these  thin 
composite  beams  (Rosakis  and  Ravi-Chandar.  1986;  Sun,  2008;  Qian  and  Sun,  2008;  Sun  and 
Qian,  2009).  On  the  other  hand,  most  of  these  beam  bending  methods  are  unable  to  eliminate 
friction  between  the  cracked  faces,  behind  the  crack  tip.  This  can  cause  a  significant  increase  in 
values  of  Kuc  which  become  higher  than  the  intrinsic  values.  For  other  non-beam  bending 
methods,  it  has  been  shown  that  the  compressive  stresses  ahead  of  the  crack  tip  significantly 
affect  the  mode  II  fracture  toughness  of  composite  materials  (Bing  and  Sun,  2007).  Also,  the 
values  of  the  measured  mode  II  fracture  toughnesses  are  strongly  dependent  on  the  type  of  test 
used  and  show  a  wide  variation.  With  the  extensive  applications  of  composite  materials,  it  is 
very  important  to  develop  new  approaches  to  measure  the  mode  II  fracture  toughnesses. 

The  losipescu  shear  test  has  been  developed  and  extensively  modified  to  measure  the  shear 
strength  of  composites  and  bonded  polymers  (Walrath  and  Adams,  1983).  This  fixture  has  been 
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used  to  characterize  the  shear  strength  of  short-beam  shear  specimens;  and  the  average  shear 
strengths  are  found  to  be  very  close  to  the  shear  strengths  of  the  losipescu  shear  specimens  in 
spite  of  very  different  interfacial  shear  stress  distributions  (Krishnan  and  Xu.  2010).  In  this 
investigation,  the  advantage  of  having  zero  interfacial  normal  stress  (zero  bending  moment)  at 
the  specimen  center  is  used  to  develop  an  efficient  experimental  method  in  determining  the  mode 
II  fracture  toughness  of  bonded  polymers.  This  method,  which  is  based  on  the  losipescu  fixture, 
minimizes  the  effect  of  friction  on  the  mode-II  fracture  toughness.  A  theoretical  formula  is 
derived  and  is  verified  using  finite  element  analysis.  This  formula  is  extended  to  include 
composite  materials  and  a  calibration  chart  is  provided.  Experimental  specimens  of  bonded 
polymers  are  tested  to  measure  fracture  toughness  values  and  the  load-displacement  curves  are 
verified  numerically  using  cohesive  element  simulation.  Indeed,  Bansal  and  Kumosa  (1998)  have 
proposed  a  “Double  Edge-Crack”  losipescu  shear  specimen  to  measure  the  mixed-mode  fracture 
toughness  of  composite  materials  subjected  to  biaxial  loads.  The  current  objective,  however,  is  to 
measure  the  mode-II  fracture  toughness  alone  with  our  approach,  which  is  much  easier  to  use 
than  previous  ones. 

2.  Theoretical  Background 

2.  1  Formula  for  the  homogenous  materials 

The  stress  intensity  factors  of  an  asymmetric  four-point  bend  specimen  (Fig.  1)  have  been 
calculated  by  Suresh  et  al.,  (1990) 
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Where  Ki  and  Kn  are  the  mode  I  and  mode  II  stress  intensity  factors  respectively;  x  is  the  shear 
stress  across  the  interface,  W  is  the  specimen  width,  a  is  the  crack  length  and  S  is  the  off-set 
distance  between  the  crack  plane  and  the  loading  point.  F|  and  Fn  are  dimensionless  functions  of 
a/W  and  are  provided  in  Suresh  et  al.  (1990).  Additionally,  the  shear  stress  per  unit  thickness  can 
be  expressed  as 

P  A-B 

T~WA  +  B  (3) 

Here,  P  is  the  applied  load,  A  and  B  are  distances  of  the  loading  points  from  the  center  as  shown 
in  Fig.l . 
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Fig.  1  Schematic  diagram  of  the  short-beam  shear  test  with  equivalent  applied 
loading  at  distances  A  and  B 

In  this  investigation,  a  short-beam  shear  specimen  with  a  bonded  interface  is  chosen  to  match 
those  dimensions  required  for  a  widely  used  Iosipescu  test  fixture  (Krishnan  and  Xu.  2010).  The 
short-beam  shear  test  is  indeed  an  anti-symmetric  four  point  bend  test.  We  choose  specific 
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adhesives  with  similar  Young's  modulus  as  the  bonded  polymers  (E=l-5GPa).  Also,  the  actual 
thickness  of  the  adhesive  interface  is  only  of  the  order  of  the  micrometers,  which  is  negligible  in 
comparison  with  other  dimensions  of  the  specimen.  Thus,  the  same-material  bonded  specimen  is 
treated  as  homogenous  materials,  however  is  not  isotropic  as  the  strength  along  the  adhesive 
interface  will  be  different  from  the  strength  of  the  bulk  material.  The  loads  on  the  Iosipescu 
fixture  are  applied  in  the  form  of  displacements  on  anti-symmetric  loading  blocks.  The 
corresponding  loads  are  not  point  loads  and  are  applied  from  the  loading  blocks  in  the  form  of 
highly  non-linear  line  loads.  However,  these  loads  can  be  expressed  as  equivalent  point  loads 
acting  at  a  specific  distance  from  the  center  line  as  shown  in  Fig.  1.  For  a  short-beam  shear 
specimen,  A  and  B  are  calculated  to  be  16  mm  and  34.5  mm  for  isotropic  materials.  These  values 
are  calculated  from  finite  element  simulations  with  the  corresponding  procedure  outlined  in 
Section  4.  In  every  case  the  moment  about  the  center  line  (along  the  interface)  is  calculated  to  be 
zero  which  is  critical  in  ensuring  that  S  from  equation  (1)  is  zero.  The  short-beam  set-up  works 
perfectly  well  for  the  case  of  the  pure  mode  11  determination  of  same-material  specimens  due  to 
the  fact  that  the  moment  about  its  center  line  is  always  zero  based  on  static  balance  analysis. 

In  the  present  case,  there  is  no  off-set  of  the  loads  from  the  center;  therefore  S  is  0.  Hence,  Ki  is  0 
and  wc  obtain  a  pure  mode  II  condition.  Further,  from  equations  (2)  and  (3),  we  obtain  a  formula 
for  Kn  as 


Kji  = 


P  A-B 
Wl  A  +  B 


(4) 


Here,  t  is  the  thickness  of  the  specimen.  Also,  in  the  present  case,  we  use  affV  of  0.5  for  the 
convenient  purpose  of  specimen  preparations.  This  gives  a  Fh  value  of  1.3  (Suresh  et  al.  1990). 
Hence,  equation  (4)  can  be  further  expressed  as 
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Kn  =  1 .3 


P  A-B 
Wt  A  +  B 


(5) 


Now,  the  energy  release  rate  Gu  can  be  expressed  as 


Gn  = 


Kn 2 

f: 


(6) 


Here,  E  is  the  effective  plane  elastic  modulus  for  both  plane  stress  and  strain  cases.  Combining 
equations  (5)  and  (6),  the  critical  energy  release  rate  Gnc  can  be  expressed  as 


Gnc  = 


1.69  Pc  2  .A-B.  2  na 

irV  ^+5*  ir 


(7) 


where  Pc  is  the  critical  load  at  the  crack  initiation  in  the  specimen.  Applying  the  values  of  A  and 
B  for  our  new  type  of  experiments  as  shown  in  Fig.l ,  we  get. 


Gnc  = 


0.7125  Pc  2  a 
E'  W  I2 


(8) 


The  applicability  of  equation  (8)  is  tested  over  a  wide  range  of  elastic  moduli  and  applied 
displacements.  This  validity  is  checked  by  obtaining  values  of  A  and  B  over  a  variety  of  elastic 
moduli  and  applied  displacements.  As  show  in  Fig.2,  values  of  A  and  B  remain  constant  for  a 
wide  range  of  Young's  moduli  (at  a  constant  applied  displacement  of  0.2  mm),  thereby  showing 
that  equation  (5)  is  valid  for  most  engineering  materials.  Fig.  3  shows  there  is  no  variation  in  A 
and  B  with  different  applied  displacement  at  a  constant  elastic  modulus  of  2.4  GPa  (Young's 
modulus  of  polymer  materials  in  this  investigation). 
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Fig.  2  Finite  element  analysis  results  showing  very  small  variations  of  A  and  B  with 
a  wide  range  of  elastic  moduli  under  an  applied  displacement  of  0.2  mm 
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Fig.  3  Finite  element  analysis  results  showing  very  small  variations  of  A  and  B  with  a  wide 
range  of  applied  displacements  for  poly  meric  specimen  with  an  elastic  modulus  of  2.4  GPa 


2.2  Formula  for  orthotropic  fibrous  composite  materials 

The  above  discussion  can  be  extended  to  orthotropic  materials  such  as  unidirectional  fiber 
composite  materials  as  shown  in  Fig.  4.  A  similar  analysis  to  determine  the  fracture  toughness  of 
composite  materials  is  conducted.  First,  the  effective  plane  elastic  modulus  of  the  composite 
materials  is  calculated  by  using  the  formula  presented  below  (Sih  et  al.,  1965:  Xu  et  al..  1996) 


1  -  1  1 

[fe  + 

2v'12  ] 

F 

^22 

e"  4ie22  ]j 

<G12 

E 

) 

2 

(9) 
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Here.  En  and  E22  are  the  elastic  moduli  of  the  composite  in  the  fiber  (longitudinal)  and 
transverse  directions.  In  this  investigation,  the  mode-II  crack  path  should  be  the  same  direction 
with  the  fiber,  or  equivalently  the  fibers  should  be  along  the  y  direction.  The  local  1-2  coordinate 
system  is  shown  in  Fig.  4.  Gn  is  the  shear  modulus  of  the  composite  and  V12  is  the  Poisson’s 
ratio  of  the  composite.  The  fracture  toughness  can  then  be  calculated  by  using  the  same 
relationship  in  equation  (6).  but  by  using  the  effective  Young's  modulus  of  the  composites.  The 
values  of  A  and  B  in  case  of  composites  are  calculated  from  finite  elements  to  be  15.5  mm  and 
32.5  mm  and  are  found  to  be  slightly  different  from  same-material  specimens.  In  order  to 
determine  values  of  A  and  B  the  loads  on  the  specimen  are  calculated  using  finite  elements.  A 
similar  approach  as  in  case  of  same-material  specimens  is  used.  The  equivalent  loading  points 
for  A  and  B  are  found  to  be  constant  and  are  plotted  against  the  ratio  of  the  typical  composite 
elastic  moduli  En  and  E22  in  Fig.  5. 
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Fig.  4  Schematic  diagram  of  a  short-beam  shear  specimen  for  unidirectional 
composites  materials 
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Fig.  5  Finite  element  analysis  results  showing  very  small  variation  in  A  and  B  with  elastic 
modulus  ratio  for  different  types  of  orthotropic  composite  materials. 


3.  Experimental  investigation 

3.1  Specimen  preparation  and  test  procedure 

Specimens  were  made  of  polycarbonate  and  PMMA.  Each  specimen  was  made  of  two  individual 
halves  of  length  38.1  mm  and  width  (W)  19.1  mm  and  was  bonded  at  the  interface  using  an 
adhesive  to  enable  interfacial  failure.  Each  specimen  had  a  thickness  of  5.5  mm.  The  bonding 
surfaces  of  the  individual  halves  were  sand-blasted  to  improve  the  bonding  quality.  While 
Weldon- 10  was  used  to  provide  strong  bonding.  Loctite  384  was  used  to  provide  weak  bonding. 
These  adhesives  were  specifically  chosen  such  that  their  elastic  modulus  when  cured  was  close 
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to  that  of  polycarbonate  and  PMMA  (2-4  GPa).  This  was  done  to  simplify  the  mechanics  of  the 
polymer/adhesive/polymer  interface.  Each  specimen  had  an  initial  crack  (a/W=0.5)  which  was 
made  by  covering  one  half  of  the  adhesive  using  a  thin  tape.  The  adhesive  was  then  applied  to 
the  other  half  and  the  specimen  was  bonded  using  a  special  fixture  to  guarantee  dimensionality. 
The  specimens  were  left  to  cure  for  a  period  of  24  hours  to  achieve  the  bonding  strength. 

The  experimental  set-up  consisted  of  three  parts  including  a  mechanical  system  to  load  the 
specimens,  an  optical  system  to  capture  fringe  patterns  and  an  imaging  set-up  to  record  the 
images  (Xu  et  al.,  2004a).  The  fringes  were  solely  used  to  observe  the  failure  process  and  to 
determine  the  crack  path.  The  mechanical  system  consisted  of  a  MTS  810  test  machine  and  an 
losipescu  fixture.  Load  was  applied  to  the  movable  part  of  the  losipescu  fixture  in  the  form  of 
displacements  at  a  rate  of  1  mm/min  until  failure.  About  seven  specimens  were  tested  to  ensure 
repeatability. 


3.2  Experimental  Results  for  Gnc 

The  experimental  results  are  presented  in  Table  1.  Two  fracture  modes  were  observed:  1)  if  the 
interfacial  bonding  was  strong,  the  initial  crack,  although  loaded  in  shear,  kinks  from  the  original 
crack  path  and  forms  a  mode-I  crack  (its  symmetrical  stress  field  was  verified  by  optical 
techniques);  and  2)  a  pure  mode-II  crack  propagates  along  the  interfacial  bonding  (a  self-similar 
crack)  if  the  interfacial  bonding  was  weak.  However,  in  case  of  the  strongly  bonded 
polycarbonate  specimens,  kinking  starts  away  from  the  crack  tip.  The  crack  was  observed 
(during  the  experiment  by  optical  methods  and  after  failure  by  observation  of  the  failure 
surfaces)  to  originate  as  a  mode  II  crack  and  then  kink  away  as  a  mode  I  crack.  Hence,  it  is 
concluded  that  the  intrinsic  mode  II  fracture  toughness  for  polycarbonate  specimens  is  obtained. 
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In  case  of  PMMA,  the  crack  was  observed  to  kink  from  the  crack  tip  itself.  So,  a  lower  bound  for 
the  Knc  value  is  obtained,  and  not  the  mode  II  fracture  toughness,  as  indicated  in  Table  1. 


Table.  1  Experimental  results  for  same-material  specimens 


Bonded 

materials 

Bonding 

Crack  initiation 
load  (N) 

K11C  (MPa  m0  5) 

Gnc(J/m2) 

Polycarbonate 

Weak 

4 1 7  ±  49 

0.33  ±  0.04 

45.8 

Polycarbonate 

Strong 

1002  ±57 

0.76  ±0.05 

260.8 

PMMA 

Weak 

451  ±  100 

0.35  ±  0.08 

34.6 

PMMA 

Strong 

975  ±81 

>0.77  ±  0.06 

>155.6 

The  values  of  Knc  are  seen  to  be  equal  or  less  than  the  values  of  Kic  for  the  same  material  types 
obtained  from  our  previous  experimental  research  (Krishnan  and  Xu,  2009).  This  indicates  that 
there  is  less  friction  betw  een  the  cracked  faces  obtained  in  our  experiments. 


3.3  Analysis  on  the  crack  kinking 

An  initial  crack  seeks  the  path  of  least  resistance  to  propagate  (Anderson,  2004).  In  the  present 
case,  a  pure  mode  II  crack  is  needed,  rather  than  a  kinked  mode  I  crack.  The  angle  of  the  crack 
kinking  can  be  calculated  theoretically  (Erdogan  and  Sih,  1963;  Williams  and  Ewing,  1972).  The 
kinking  angle  has  been  calculated  for  a  mixed  mode  case  by  Bhattacharjee  and  Knott  (1995).  A 
schematic  diagram  of  crack  kinking  along  with  the  conventions  is  shown  in  Fig.  6.  The  criterion 
for  an  initial  crack  to  kink  along  an  angle  can  be  expressed  as 
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Where  the  shear  stress  is  zero  and  the  hoop  stress  is  maximized.  The  hoop  stress  and  shear  stress 
close  to  a  crack  tip  can  be  expressed  in  polar  coordinates  as  (Anderson,  2004) 
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The  condition  in  equation  (10)  is  solved  for  0  and  the  predicted  angle  of  the  kink  is  calculated  as 
19.5°  from  the  x-axis.  This  is  close  to  the  angle  of  kink  obtained  experimentally  from 
polycarbonate  and  PMMA.  It  should  be  noted  that  the  kinking  angle  is  independent  of  the 
material  used  if  the  T-stress  is  neglected.  The  effect  of  the  T-stress  before  and  after  crack  kinking 
has  been  reported  by  Li  and  Xu  (2007).  A  typical  strongly  bonded  polycarbonate  specimen 
showing  a  kinking  angle  of  about  25°  from  the  x-axis  is  presented  in  Fig.  7. 

In  order  to  predict  which  failure  mode  will  occur  at  first,  the  energy  release  rates  and  fracture 
toughnesses  of  two  failure  modes  should  be  compared  (Hutchinson  and  Suo,  1992;  Xu  and 
Rosakis,  2003;  Xu  et  ah,  2003).  The  energy  release  rate  of  the  initial  crack  is  a  function  of  the 
potential  crack  path  or  a  function  of  0.  For  the  kinked  mode-I  crack,  the  energy  release  rate, 

Y  pm 

G(0=0C)  should  exceed  the  mode-I  fracture  toughness  of  the  polymer  material  1 IC  ,  i.e., 

Gi0  =  0cKl  (13) 

yPM  - 

1  Ic 
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Where  the  crack  kinking  angle  0C  is  around  20°  from  the  x  axis  based  on  Erdogan  and  Sih  (1963). 
For  the  pure  mode-11  crack  initiation  and  propagation,  the  energy  release  rate,  G(0=9O°)  should 

F/r  . 

exceed  the  mode-11  fracture  toughness  of  the  interface  llc ,  i.e., 

G(6>  =  90°)  (14) 

rIT 

1  Ik 

So,  the  competition  of  the  energy  release  rate  and  the  fracture  toughness  leads  to  different  failure 
modes.  Here,  if  equation  (13)  is  satisfied  at  first,  a  crack  loaded  in  shear  will  kink  as  a  mode-1 
crack.  If  equation  (14)  is  satisfied  at  first,  a  crack  loaded  in  shear  will  propagate  as  a  mode-11 
crack  along  the  interface. 
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Fig.  7  Picture  of  a  typical  strongly  bonded  polycarbonate  specimen 
showing  crack  kinking  from  the  interface 

4.  Numerical  Methods 

4.1  Finite  element  model  to  compare  with  formula  for  same-material  specimens 

A  two-dimensional  Finite  element  analysis  was  conducted  using  Abaqus®  6.9.1.  The  mesh  is 
progressively  graded  in  order  to  reduce  computational  cost,  and  is  refined  near  the  crack  tip  to 
best  capture  the  stress  singularity.  Polycarbonate  specimens  with  an  elastic  modulus  of  2.4  GPa 
and  a  Poisson’s  ratio  of  0.37  are  chosen  for  the  analysis.  The  dimensions  of  the  model  are  the 
same  as  the  experimental  specimens  as  described  in  Section  2.  The  geometry  of  each  individual 
part  is  symmetric  while  the  loading  on  the  material  is  anti-symmetric.  The  model  is  assumed  to 
be  monolithic  with  the  zero  thickness  of  the  adhesive  layer. 

The  loading  is  applied  in  the  form  of  loading  block  on  the  specimen  edges  similar  to  what  is 
observed  in  reality.  In  order  to  incorporate  a  realistic  simulation  of  the  loads  on  the  specimen  an 
iterative  procedure  is  adopted  (Xu  et  al.,  2004b;  Krishnan  and  Xu.  2010).  The  loads  in  the  form 
of  displacements  are  applied  to  the  movable  part  of  the  specimen  while  the  other  part  is  held 
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fixed.  After  each  analysis  the  reaction  forces  on  the  loading  edges  are  checked  to  see  if  they  are 
in  compression  or  not,  because  in  experiments,  only  compressive  loading  is  observed.  The 
constraints  are  removed  from  those  nodes  which  showed  tensile  reaction  forces  and  the  analysis 
is  repeated  until  convergence.  The  stress  intensity  factors.  K]  and  Kn  are  calculated  from  the 
finite  element  analysis  and  are  compared  with  the  results  from  equation  (5).  The  results  are 
presented  in  Table  2.  It  can  be  seen  that  the  Ki  is  negligibly  small  in  comparison  with  Kn  as  is 
expected  from  the  theoretical  results  presented  earlier.  Also,  the  values  of  Kn  from  finite  element 
analysis  agree  well  with  the  values  calculated  from  our  formula  for  same-material  specimens. 
This  validates  our  theoretical  formula  for  same-material  specimens. 


Table  2.  Variation  of  mode  I  and  mode  I  stress  intensity  factors  Ki  and  Kn  with  applied 
load  for  same-material  polycarbonate  specimen 


Load  (N) 

Finite  Element  Analysis 

Formula 

K,  (MPa  m0-5) 

K„  (MPa  m0-5) 

K,  (MPa  m0  5) 

Kn  (MPa  m0-5) 

135 

0.0019 

0.1067 

0 

0.1060 

405 

0.0057 

0.3199 

0 

0.3179 

675 

0.0097 

0.5333 

0 

0.5299 

1350 

0.019 

1.067 

0 

1.0598 

A  snapshot  of  the  deformed  mesh  near  the  crack  tip  at  failure  load  is  presented  in  Fig.  8.  It  can  be 
seen  that  there  is  no  crack  surface  separation  (opening  mode-I),  and  only  sliding  along  the  crack 
face  (shear  mode-II).  This  confirms  that  there  is  no  mode  I  component  from  our  finite  element 
analysis.  The  shear  and  normal  stresses  along  the  interface  at  a  failure  load  of  1 000  N  are  plotted 
in  Fig.  9.  The  graph  demonstrates  the  singularity  of  the  stresses  at  the  crack  tip.  Also,  there  is  no 
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contact  between  the  crack  lips  indicating  that  there  is  little  friction  involved  during  experiments. 
Fig.  10  shows  a  comparison  of  shear  stresses  along  the  crack  faces  for  different  loading  cases. 
These  values  are  seen  to  be  close  to  zero  even  at  higher  loads  (beyond  failure  loads) 
demonstrating  the  absence  of  friction  between  the  crack  faces. 


Shear  Deformation 

Fig.8  Snapshot  of  deformed  mesh  (near  the  crack  tip)  at 
failure  load.  This  mesh  is  used  for  modeling  short-beam 
shear  specimen 
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Fig.9  Variation  of  shear  stress  and  normal  stress  along  the  interface  for  polycarbonate  same- 
material  bonding  at  a  failure  load  of  1 000  N. 


18 


Krishnan  and  Xu  (2010) 


Efficient  approach  to  measure  mode-II  fracture  toughness 


Shear  Stress  (MPa) 


Fig.  10  Variation  of  shear  stresses  along  the  interface  for  polycarbonate  same-material  bonding 
with  increase  in  applied  loading.  There  is  negligible  friction  observed  beyond  the  crack  tip. 


4.2  Numerical  modeling  for  composites 

A  similar  model  as  described  in  above  is  made  for  composite  materials  with  Fiber  direction  along 
the  y-direction  as  shown  in  Fig.4.  This  type  of  unidirectional  composite  materials  ensures  that 
the  mode-II  crack  will  not  kink  and  will  propagate  along  the  Fiber-matrix  interface  due  to  strong 
constraint  of  the  fiber.  Three  kinds  of  typical  composite  material  data  are  chosen  from  Daniel 
and  Ishai  (2005)  to  conduct  the  numerical  analysis  (Table  3). 

However,  in  case  of  orthotropic  unidirectional  composite  materials,  we  do  not  have  the  graph  of 
Fn  which  has  been  provided  for  same-material  specimens  by  Suresh  et  al.  (1990)  and  for  bi- 
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material  specimens  by  O'Dowd  et  al.  (1992).  Hence,  the  plot  of  P'n  as  a  function  of  a/W  is 
calculated  for  three  different  types  of  composite  materials,  viz.  E-Glass/Epoxy  (55  %  volume 
percent  of  fibers,  unidirectional  GFRP),  Kevlar  49/Epoxy  (60%  volume  percent  of  fibers, 
unidirectional  KFRP)  and  AS4/3501-6  Carbon/Epoxy  (63%  volume  percent  of  fibers, 
unidirectional  CFRP).  The  plot  of  Fn  as  a  function  of  a/W  is  presented  in  Fig.  1 1.  It  is  seen  that 
this  plot  is  almost  constant  for  these  three  types  of  typical  composites.  While  calculating  the 
mode  II  stress  intensity  factors  for  composites,  this  chart  becomes  necessary  to  use  the  current 
approach  and  is  complementary  to  the  one  presented  in  Suresh  et  al.  (1990).  Also,  the 
numerically  calculated  values  of  Ki,  Kn  and  Fn  at  different  values  of  a/W  ratios  are  presented  in 
Table  4  for  the  three  typical  composite  types  chosen.  The  ratio  of  Ki/Kn  remains  negligibly  small 
indicating  that  the  crack  is  indeed  a  pure  mode  II  crack  in  unidirectional  composites.  The  applied 
loads  on  each  composite  specimen  (corresponding  to  a  fixed  applied  displacement  of  0.5  mm) 
are  also  presented  in  Table  4.  The  different  types  of  composites  take  different  loads  (in  spite  of  a 
fixed  applied  displacement)  due  to  their  different  elastic  properties.  However,  the  dimensionless 
stress  intensity  factor,  Fn,  remains  very  similar  for  the  different  kinds  of  composites. 

Table  3.  Properties  for  three  typical  composites  (Daniel  and  Ishai,  2005) 


Materials 

Fiber  volume 
percent  (%) 

Eji 

(GPa) 

E22 

(GPa) 

G12 

(GPa) 

G23 

(GPa) 

V]2 

V23 

Carbon/Epoxy 

63 

147 

10.3 

7 

3.7 

0.27 

0.54 

E-Glass/Epoxy 

55 

41 

10.4 

4.3 

3.5 

0.28 

0.5 

Kevlar/Epoxy 

60 

80 

5.5 

2.2 

1.8 

0.34 

0.4 
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Table  4.  Variation  of  Fn  w  ith  a/W  for  different  types  of  composites.  Ki  and  Kh  are  expressed 
in  MPa  m]  7  (A=l  5.49  mm,  B=32.44  mm) 


CFRP 


a/W 

Load  (N) 

Ki  (MPa  m1/2) 

Kh  (MPa  m1/2) 

Kn/Ki 

Fn 

0.2 

1009 

0.0124 

1.2412 

0.0099 

0.3836 

0.3 

1008 

0.0090 

2.1345 

0.0042 

0.6603 

0.4 

1004 

0.0026 

3.0042 

0.0009 

0.933 

0.5 

1001 

0.0006 

3.8580 

0.0001 

1.2018 

0.6 

991 

0.0083 

4.7118 

0.0018 

1.4825 

0.7 

979 

0.0146 

5.6605 

0.0026 

1.8029 

GFRP 


a/W 

Load  (N) 

K]  (MPa  m,/J) 

Kn  (MPa  ml/J) 

K„/K, 

F„ 

0.2 

709 

0.0405 

0.9645 

0.0419 

0.4242 

0.3 

708 

0.0281 

1 .6286 

0.0173 

0.7173 

0.4 

704 

0.0072 

2.2927 

0.0031 

1.0155 

0.5 

698 

0.0095 

2.9410 

0.0032 

1.3138 

0.6 

687 

0.0263 

3.5418 

0.0075 

1.6075 

0.7 

671 

0.0361 

4.2058 

0.0086 

1.9544 

KFRP 


a/W 

Load  (N) 

Ki  (MPa  m1/2) 

Kn  (MPa  m1/J) 

K„/K, 

Fn 

0.2 

420 

0.0058 

0.5139 

0.0113 

0.3815 

0.3 

420 

0.0039 

0.8854 

0.0044 

0.6573 

0.4 

418 

0.0006 

1.2491 

0.0005 

0.9318 

0.5 

417 

0.0004 

1.6128 

0.0003 

1.2059 

0.6 

413 

0.0043 

1.9606 

0.0022 

1.4802 

0.7 

408 

0.0063 

2.3401 

0.0027 

1.7884 
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Fig.  1 1  Variation  of  dimensionless  stress  intensity  factor  Fn  with  relative  crack  depth  for  different 
ty  pes  of  composites. 


4.3  Cohesive  element  simulation 

In  order  to  simulate  the  experimental  load-displacement  curves  and  the  failure  load,  a  cohesive 
element  analysis  is  employed.  Cohesive  elements  are  exclusively  used  along  the  bonded 
interface  in  order  to  model  the  failure  process.  Cohesive  elements  with  zero  thickness  are  used  in 
the  model.  A  bi-linear  cohesive  law  is  considered  for  this  case  and  a  representative  law  for 
strongly  bonded  polycarbonate  is  shown  in  Fig.  12.  Such  types  of  cohesive  laws  have  been  used 
successfully  in  the  past  to  model  mixed  mode  failure  in  brittle  materials  (Camacho  and  Ortiz, 
1996;  Ruiz  et  al..  2001).  The  four  types  of  specimens  in  Table  1  are  modeled  using  a  similar 
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bilinear  cohesive  law  but  with  different  values  of  constants  (shown  in  Fig.  12).  The  value  of 
fracture  toughness,  Ghc,  obtained  from  the  current  experimental  results  (Table  1)  is  used  here 
and  equals  the  total  area  enclosed  by  the  cohesive  law.  The  elastic  modulus  of  the  bulk  material 
is  used  as  the  slope  of  the  initial  linear  part  of  the  traction  separation  law  .  The  maximum  value  of 
traction  (t)  is  obtained  from  our  previous  measurements  for  interfacial  shear  strengths  of  bonded 
polycarbonate  systems  (Krishnan  and  Xu,  2009).  The  critical  opening  displacement  (5C)  can  then 
be  calculated  from  these  known  values.  The  finite  element  model  used  is  similar  to  the  one 
described  in  the  previous  sections  except  in  the  inclusion  of  the  cohesive  elements.  Loads  are 
then  applied  to  the  finite  element  model  as  described  previously.  The  simulated  load- 
displacement  curves  are  plotted  along  with  the  experimental  load-displacement  curves  in  Fig.  13. 
Three  different  experimental  curves  are  shown  to  depict  the  variation  in  the  experiments.  The 
cohesive  element  model  is  seen  to  capture  the  initial  slope  of  the  experimental  load-displacement 
curves.  There  is  only  a  reasonable  amount  of  error  for  the  value  of  the  final  failure  load  as  shown 
by  the  cohesive  model.  A  similar  analysis,  albeit  with  a  different  cohesive  law,  is  carried  out  for 
weakly  bonded  polycarbonate  specimens  and  the  corresponding  results  are  shown  in  Fig.  14. 
Different  values  of  fracture  toughness  and  maximum  traction  are  used  in  this  case  and  these  data 
were  from  our  previous  measurement  for  the  same  material  system.  This  model  predicts  the 
crack  initiation  load  and  failure  displacement  with  a  reasonable  amount  of  error.  The  graphs  in 
Figs.  13  and  14  are  plotted  on  the  same  scale  to  allow  for  a  fair  comparison  between  the  crack 
initiation  loading  of  the  strong  and  the  weak  bonds. 
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Displacement  (mm) 


Fig.  12  Traction-Separation  law  used  for  strongly  bonded  Polycarbonate 
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Fig.  13  Simulated  and  experimental  load-displacement  graphs  for  strongly  bonded  polycarbonate 
systems. 


The  load-displacement  curves  obtained  from  the  proposed  new  approach  show  a  sudden  and 
clear  drop  as  seen  in  Figs.  13  and  14.  The  crack  initiation  and  the  final  failure  loads  coincide. 
This  is  very  accurate  in  finding  the  crack  initiation  load  and  is  a  major  advantage  of  our  new 
approach.  For  other  methods  involving  beam  bending,  the  crack  initiation  load  is  not  obvious  as 
the  beam  compliance  change  involving  short  crack  propagation  is  small  (Xu  and  Kou,  1994). 
This  is  a  possible  reason  for  overestimation  of  the  mode  II  fracture  toughness  of  composite 
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materials.  Future  work  involves  extending  this  approach  to  measure  the  fracture  toughness  of 
unidirectional  composites. 


Fig.  14  Simulated  and  experimental  load-displacement  graphs  for  weakly  bonded 
polycarbonate  systems. 


5.  Conclusions 

The  mode  II  fracture  toughness  of  materials  with  preferred  interfaces  can  be  measured  by  using 
the  proposed  short-beam  shear  approach.  The  amount  of  friction  between  the  cracked  surfaces 
behind  the  crack  tip  is  negligible,  and  hence  an  intrinsic  value  of  the  pure  mode-1 1  fracture 
toughness  is  obtained.  Another  feature  of  our  new  approach  is  the  accurate  measurement  of  crack 
initiation  load  which  is  not  available  in  previous  techniques.  The  current  method  is  justified  by 
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the  use  of  numerical  techniques  (finite  element  analysis  and  cohesive  element  analysis),  and 
comparison  with  experimental  data.  This  method  can  also  be  extended  to  unidirectional 
composite  materials. 
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Fig.  4  Schematic  diagram  of  a  short-beam  shear  specimen  for  unidirectional  composites 
materials 

Fig.  5  Finite  element  analysis  results  showing  very  small  variation  in  A  and  B  w  ith  elastic 
modulus  ratio  for  different  types  of  orthotropic  composite  materials. 

Fig.  6  Schematic  diagram  showing  the  crack  kinking 

Fig.  7  Picture  of  a  ty  pical  strongly  bonded  polycarbonate  specimen  showing  crack  kinking  from 
the  interface 

Fig. 8  Snapshot  of  deformed  mesh  (near  the  crack  tip)  at  failure  load.  This  mesh  is  used  for 
modeling  short-beam  shear  specimen 

Fig. 9  Variation  of  shear  stress  and  normal  stress  along  the  interface  for  polycarbonate  same- 
material  bonding  at  a  failure  load  of  1000  N. 

Fig.  10  Variation  of  shear  stresses  along  the  interface  for  polycarbonate  same-material  bonding 
with  increase  in  applied  loading.  There  is  negligible  friction  observed  beyond  the  crack  tip. 

Fig.  11  Variation  of  dimensionless  stress  intensity  factor  Fn  w  ith  relative  crack  depth  for 
different  types  of  composites. 

Fig.  12  Traction-Separation  law  used  for  strongly  bonded  Polycarbonate 

Fig.  13  Simulated  and  experimental  load-displacement  graphs  for  strongly  bonded  polycarbonate 
systems. 

Fig.  14  Simulated  and  experimental  load-displacement  graphs  for  w  eakly  bonded  polycarbonate 
systems. 
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Fig.  1  Schematic  diagram  of  the  short-beam  shear  test  with  equivalent  applied 
loading  at  distances  A  and  B 
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Fig.  2  Finite  element  analysis  results  showing  very  small  variations  of  A  and  B  with 
a  w  ide  range  of  elastic  moduli  under  an  applied  displacement  of  0.2  mm 
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Fig.  3  Finite  element  analysis  results  showing  very  small  variations  of  A  and  B  with  a  wide 
range  of  applied  displacements  for  polymeric  specimen  with  an  elastic  modulus  of  2.4  GPa 
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Fig.  5  Finite  element  analysis  results  showing  very  small  variation  in  A  and  B  with  elastic 
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Fig.  7  Picture  of  a  typical  strongly  bonded  polycarbonate  specimen 
showing  crack  kinking  from  the  interface 
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Shear  Deformation 

Fig.8  Snapshot  of  deformed  mesh  (near  the  crack  tip)  at 
failure  load.  This  mesh  is  used  for  modeling  short-beam 
shear  specimen 
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Fig.9  Variation  of  shear  stress  and  normal  stress  along  the  interface  for  polycarbonate  same- 
material  bonding  at  a  failure  load  of  1000  N. 
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Fig.  10  Variation  of  shear  stresses  along  the  interface  for  polycarbonate  same-material  bonding 
with  increase  in  applied  loading.  There  is  negligible  friction  observed  beyond  the  crack  tip. 
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Fig.  11  Variation  of  dimensionless  stress  intensity  factor  Fn  with  relative  crack  depth  for  different 
types  of  composites. 
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Fig.  12  Traction-Separation  law  used  for  strongly  bonded  Polycarbonate 
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Fig.  13  Simulated  and  experimental  load-displacement  graphs  for  strongly  bonded  polycarbonate 
systems. 
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Fig.  14  Simulated  and  experimental  load-displacement  graphs  for  weakly  bonded  polycarbonate 
systems. 
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Table.  1  Experimental  results  for  same-material  specimens 


Bonded 

materials 

Bonding 

Crack  initiation 
load  (N) 

K„c  (MPa  m0'5) 

G„c(J/m2) 

Polycarbonate 

Weak 

41 7  ±49 

0.33  ±  0.04 

45.8 

Polycarbonate 

Strong 

1002  ±57 

0.76  ±0.05 

260.8 

PMMA 

Weak 

451  ±  100 

0.35  ±0.08 

34.6 

PMMA 

Strong 

975  ±81 

>0.77  ±  0.06 

>155.6 
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Table  2.  Variation  of  mode  I  and  mode  I  stress  intensity  factors  Ki  and  Kn  with  applied 
load  for  same-material  polycarbonate  specimen 


Load  (N) 

Finite  Element  Analysis 

Formula 

Ki  (MPa  m°  5) 

Kn  (MPa  m°  5) 

Ki  (MPa  m0  5) 

K„  (MPa  m0-5) 

135 

0.0019 

0.1067 

0 

0.1060 

405 

0.0057 

0.3199 

0 

0.3179 

675 

0.0097 

0.5333 

0 

0.5299 

1350 

0.019 

1.067 

0 

1.0598 
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Table  3.  Properties  for  three  typical  composites  (Daniel  and  Ishai,  2005) 


Materials 

Fiber  volume 
percent  (%) 

En 

(GPa) 

E  22 
(GPa) 

Gn 

(GPa) 

G23 

(GPa) 

VJ2 

V23 

Carbon/Epoxy 

63 

147 

10.3 

7 

3.7 

0.27 

0.54 

E-Glass/Epoxy 

55 

41 

10.4 

4.3 

3.5 

0.28 

0.5 

Kevlar/Epoxy 

60 

80 

5.5 

2.2 

1.8 

0.34 

0.4 
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Table  4.  Variation  of  Fn  with  a/W  for  different  types  of  composites.  Ki  and  Kn  are  expressed 
in  MPa  m1 2  (A=15.49  mm,  B=32.44  mm) 


CFRP 


a/W 

Load  (N) 

Ki  (MPa  m,/2) 

Kn  (MPa  m1/2) 

K„/K, 

F„ 

0.2 

1009 

0.0124 

1.2412 

0.0099 

0.3836 

0.3 

1008 

0.0090 

2.1345 

0.0042 

0.6603 

0.4 

1004 

0.0026 

3.0042 

0.0009 

0.933 

0.5 

1001 

0.0006 

3.8580 

0.0001 

1.2018 

0.6 

991 

0.0083 

4.7118 

0.0018 

1.4825 

0.7 

979 

0.0146 

5.6605 

0.0026 

1.8029 

GFRP 


a/W 

Load  (N) 

Ki  (MPa  m,/2) 

Kn  (MPa  m1/2) 

K„/K, 

F„ 

0.2 

709 

0.0405 

0.9645 

0.0419 

0.4242 

0.3 

708 

0.0281 

1.6286 

0.0173 

0.7173 

0.4 

704 

0.0072 

2.2927 

0.0031 

1.0155 

0.5 

698 

0.0095 

2.9410 

0.0032 

1.3138 

0.6 

687 

0.0263 

3.5418 

0.0075 

1.6075 

0.7 

671 

0.0361 

4.2058 

0.0086 

1.9544 

KFRP 


a/W 

Load  (N) 

Ki  (MPa  m1/2) 

Kn  (MPa  m1'2) 

K„/K, 

F„ 

0.2 

420 

0.0058 

0.5139 

0.0113 

0.3815 

0.3 

420 

0.0039 

0.8854 

0.0044 

0.6573 

0.4 

418 

0.0006 

1.2491 

0.0005 

0.9318 

0.5 

417 

0.0004 

1.6128 

0.0003 

1.2059 

0.6 

413 

0.0043 

1.9606 

0.0022 

1 .4802 

0.7 

408 

0.0063 

2.3401 

0.0027 

1.7884 
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Abstract 

A  systematic  experimental  investigation  to  determine  the  shear,  tensile  and  fracture  properties  of 
adhesive  joints  with  bonded  same-materials  (polymer-polymer)  and  bi-materials  (metal-polymer) 
is  reported.  Full-field  optical  techniques  including  photoelasticity  and  Coherent  Gradient  Sensing 
(CGS)  are  employed  to  record  the  stress  development  and  failure  in  these  adhesive  joints.  Five 
types  of  strong  and  weak  adhesives  are  used  in  conjunction  with  five  different  types  of  materials 
(aluminum,  steel.  PMMA,  polycarbonate  and  Homalite)  to  produce  a  variety  of  bonded  material 
systems.  Weldon- 10  and  Polyester  bonding  consistently  show  higher  tensile  and  shear  bonding 
strengths.  However,  bi-material  fracture  toughnesses  tend  to  show  a  lesser  difference  between 
the  weak  and  strong  adhesive  bonding.  Bi-material  systems  in  shear  and  fracture  report  lower 
properties  than  the  same-material  systems  due  to  a  higher  property  mismatch  in  the  former.  The 
resulting  complete  experimental  data  are  expected  to  be  immensely  helpful  to  computational 
mechanists  in  simulating  failure  mechanics  of  adhesive  joints. 
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1.  INTRODUCTION 

A  dissimilar-material  joint  is  a  special  type  of  material  boundary  at  which  two  or  more 
different  materials  are  joined  together  (e.g.,  aluminum  with  polycarbonate)  by  means  of  an 
adhesive.  This  type  of  an  interface  poses  challenging  problems  in  characterizing  its  associated 
mechanical  properties,  especially  if  the  constituent  materials  have  a  significantly  high  property 
mismatch.  Recently,  dissimilar-material  joints  have  found  several  engineering  applications. 
Examples  include  adhesively  bonded  composite  and  metal  joints  in  advanced  aircraft  and  ship 
structures,  fiber/matrix  interfaces  in  composite  materials,  thin  film/substrate  structure  in  micro¬ 
electromechanical  systems  (MEMS)  among  others.  However,  a  variety  of  studies  have  indicated 
that  failure  often  occurs  along  the  interface  between  two  different  materials  with  high  property 
mismatch  (e.g.,  free-edge  delamination  in  composite  laminates,  debonding  between  thin 
films/substrates),  and  that  improving  the  interfacial  mechanical  properties  (especially  reducing 
the  interfacial  stress  level)  can  modify  overall  material/structural  behavior  (Hutchinson  and  Suo, 
1992;  Gupta  et  al.,  1994;  Lara-Curzio  et  al.,  1995;  Needleman  and  Rosakis,1999;  Martin  et  al., 
2001;  Xu  and  Rosakis,  2002;  Wang,  2007).  However,  the  interfacial  strength  measurement  of 
dissimilar-material  joints  remains  a  challenge  due  to  the  stress  singularity  problem  (Reedy  and 
Guess,  1993;  Tandon  et  al.,  1999;  Akisanya  and  Fleck,  1997),  i.e.,  the  theoretical  interfacial 
stress  will  be  infinite  at  the  free-edges.  On  the  other  hand,  modem  numerical  tools  like  the 
cohesive  element  method  have  an  urgent  need  for  interfacial  strengths  and  fracture  toughnesses 
as  important  input  data.  Hence,  it  becomes  necessary  to  develop  reliable  quantitative 
measurements  in  order  to  characterize  mechanical  properties  of  dissimilar-material  joints. 
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The  primary  objective  of  this  paper  is  to  measure  complete  mechanical  properties  for 
polymer  and  metal  combinations  with  different  adhesive  bonding.  The  mechanics  of  bi¬ 
materials  is  complicated;  therefore  same-material  adhesive  joints  without  high  property 
mismatch  will  be  used  to  highlight  certain  simple  mechanics  aspects  involved  in  the  problem. 
Additionally,  we  utilize  new  joint  designs  to  remove  the  free-edge  stress  singularity,  thereby 
providing  reasonable  tensile  bonding  strength  measurements  of  dissimilar-material  joints.  We 
also  present  resulting  fringe  patterns  from  the  Coherent  Gradient  Sensing  (CGS)  and 
photoelasticity  experiments  (Singh  et  al.,  1997),  which  aid  in  validating  our  finite  clement 
model.  For  the  sake  of  completeness  and  direct  comparison,  we  include  some  mechanical 
property  data  from  our  previous  measurements  of  dissimilar-material  joint  systems. 


2.  EXPERIMENTAL  TECHNIQUE  AND  PROCEDURES 

2. 1  Jointed  material  combinations  w  ith  strong  and  w  eak  adhesive  bonding 

This  section  describes  the  preparation  of  specimens  and  related  experimental  procedures. 
Specimens  were  made  of  five  different  types  of  materials  including  aluminum,  steel, 
polycarbonate,  PMMA,  and  Homalite.  While  the  same-material  joint  specimens  were  made  of 
PMMA,  Homalite,  polycarbonate  and  aluminum,  bi-material  specimen  pairs  included 
PMMA/aluminum,  PMMA/steel.  polycarbonate/aluminum,  and  Homalite/steel.  Two  major  types 
of  adhesives  were  used  as  bonding  agents  in  most  of  the  specimens.  Loctite-384  (Henkel  Corp.. 
Rocky  Hill.  CT),  a  type  of  tough  Acrylic,  was  used  to  provide  weak  bonding.  Weldon- 10  (IPS 
Corp.,  Gardena,  CA),  a  mixture  of  Acrylic  Resin  and  Methyl  Methacrylate  Monomer,  was  used 
to  provide  strong  bonding.  Polyester,  Loctite  5083  and  Loctite  330  were  also  used  as  adhesives 
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for  interfacial  bonding.  All  these  adhesives  were  specifically  chosen  because  their  Young’s 
modulus  when  cured  is  close  to  that  of  PMMA,  polycarbonate  and  Homalite  (around  2-4GPa). 
This  ensures  that  the  mechanics  of  metal/adhesive/polymer  joints  remains  a  bi-material 
mechanics  problem  and  does  not  become  a  complicated  case  involving  three  different  kinds  of 
materials.  Every  specimen  was  bonded  together  from  separate  halves  to  enable  an  interfacial 
failure  and  all  of  the  individual  bonding  surfaces  were  sandblasted  prior  to  bonding  in  order  to 
improve  the  bonding  quality.  The  bonding  itself  was  enabled  by  a  special  fixture  to  guarantee 
dimensionality.  The  specimens  were  cured  for  a  period  of  24  hours  in  order  to  achieve  the 
bonding  strength. 


Fig.  1 

2. 2  Two  optical  techniques  and  experimental  procedures 

In  order  to  record  in-situ  fringe  patterns  developed  during  the  loading  process,  a 
mechanical-optical  system  as  illustrated  in  Fig.  1  was  used.  This  setup  included  an  optical 
system,  a  mechanical  testing  system  and  an  imaging  system.  The  mechanical  testing  system 
consisted  of  a  MTS  810  test  machine  and  a  fixture.  The  optical  system  was  utilized  to  capture 
the  in-situ  fringe  pattern  development  during  tests.  A  laser  beam  was  transmitted  through  the 
transparent  polymer  specimens,  and  the  resulting  fringe  patterns  were  recorded  by  a  digital 
camera.  The  isochromatic  fringe  patterns  observed  in  polycarbonate  and  Homalite  specimens  are 
the  contours  of  the  maximum  in-plane  shear  stress. 
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where  Gi  and  02  are  in-plane  principal  stresses,  N  is  the  fringe  order,  fa  is  the  stress-fringe 
constant  and  h  is  the  specimen  thickness  (Kobayashi,  1987).  The  optical  system  included  a  He- 
Ne  laser  source  (17  mW),  a  laser  collimator,  and  a  reflection  mirror.  The  collimator  was  used  to 
provide  a  large  and  collimated  laser  beam  of  approximately  50  mm  diameter.  The  purpose  of  the 
mirror  was  to  adjust  the  laser  beam  to  a  desired  position  for  a  specific  experiment.  The  imaging 
system  included  a  digital  camera  used  to  capture  fringe  development  and  a  density  filter  in  front 
of  the  camera  to  reduce  the  laser  intensity  of  the  laser  beam  entering  the  camera  directly.  A 
convex  lens  with  a  focal  length  150  mm  w'as  added  to  the  system  to  record  the  whole  field  of 
view.  An  important  issue  in  obtaining  good-quality  images  is  to  focus  the  digital  camera  at 
infinity,  and  to  ensure  the  distance  between  the  convex  lens  and  the  specimen  to  be  slightly 
larger  than  the  focal  length  of  the  convex  lens.  Another  technique  used  in  this  study.  Coherent 
Gradient  Sensing  technique  developed  by  Tippur  et  al.  (1991),  was  used  to  obtain  fringes  of 
gradients  of  01+02  in  PMMA  specimens  as  seen  in  Fig.2. 

Fig.2 


3.  TENSILE  STRENGTHS  OF  SAME  AND  DISSIMILAR-MATERIAL  JOINTS 

3. 1  Free-edge  stress  singularity  and  specimen  designs  to  remove  the  stress  singularity 

For  some  specific  bi-material  comers  or  edges,  several  researchers  including  Williams 
(1952),  Bogy  (1971),  have  shown  that  stress  singularities  exist.  The  asy  mptotic  stress  field  of  a 
bi-material  comer  can  be  expressed  by 

.v 

<r„(r,e)  =  E  r+KJ^O)  ( Lj  =  1,2,3)  (2) 

k-0 
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where  fjjk  (0)  is  an  angular  function  and  Kk  is  also  known  as  the  “stress  intensity  factor”.  The 
fracture  mechanics  terminology  “stress  intensity  factor”  is  used  in  interfacial  mechanics  to 
characterize  a  similar  stress  singularity  problem.  It  should  be  noticed  that  for  an  interfacial 
fracture  problem  (assuming  initial  debonding),  the  stress  singularity’  at  a  crack  tip  is  intrinsic  and 
cannot  be  removed.  However,  the  stress  singularity  in  interfacial  strength  investigation  such  as 
bi-material  joints  (assuming  perfect  bonding)  can  be  removed  through  appropriate  designs.  The 
stress  singularity  order  L  may  be  real  or  complex,  and  the  singularity  order  is  -0.5  for  a  crack  in 
the  same  material  based  on  the  Linear  Elastic  Fracture  Mechanics  (LEFM).  Hence  the  theoretical 
stress  values  will  become  infinite  as  r  (defined  in  Fig.  2(b))  approaches  zero,  if  X  has  a  positive 
real  part.  This  leads  to  a  problem  referred  to  as  the  “free-edge  stress  singularity  problem”.  It  is 
the  presence  of  this  stress  singularity  that  leads  to  erroneous  results  in  current  interfacial  strength 
measurements  besides  being  responsible  for  free-edge  debonding  or  delamination  in  dissimilar 
material  joints.  But.  if  k  has  a  non-positive  real  part,  then,  the  stress  singularity  disappears. 

Bogy  (1971)  found  that  the  stress  singularity  was  purely  determined  by  the  material 
property  mismatch  and  two  joint  angles  of  the  bi-material  comer  dx,02 (defined  in  Fig.  2(b)). 
Generally,  the  material  property  mismatch  can  be  expressed  in  terms  of  the  Dundurs’  parameters 
a  and  /?-  two  non-dimensional  parameters  computed  from  elastic  constants  of  two  bonded 
materials  (Hutchinson  and  Suo,  1992): 

H\m2  -/v»,  „  -2)-/ty(Wi -2) 

a  = -  B  = -  (5) 

Hxm2  fixm2+fi2mi 

where  pLx  is  the  shear  modulus  of  material  1,  ji2  is  the  shear  modulus  of  material  2, 
m  =  4(l-v)for  plane  strain,  v  is  the  Poisson’s  ratio  and  m  =  4/(1  +v)  for  generalized  plane 
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stress.  The  stress  singularity  order  is  related  to  material  and  geometric  parameters,  and  is 
determined  by  a  characteristic  equation  of  coefficients  A  (6|,  02,  p)  through  F  (0j,  O2.  p): 

f(e„82,a.p.p)  =  Ap2+2Bap  +  Ca2+2Dp  +  2Ea  +  F  =  0  (4) 

where  p=l-A  and  A,  B,  C,  D,  E  and  F  are  as  defined  by  Bogy  (1971).  Therefore,  our  basic  idea 
is  to  vary  these  four  independent  parameters (#,,#2, a,p)  in  order  to  obtain  a  negative  real  part 
of  the  stress  singularity  order/.  Xu  et  al.,  (2004a)  chose  an  interfacial  design  with  two  joint 
angles:  8i  =  65°  and  02  =  45*’  and  assumed  material  1  to  be  a  typical  hard  material  and  material 
2  to  be  a  soft  material.  Then,  there  will  be  no  stress  singularity  for  a  wide  range  of  current 
engineering  materials  (small  deviation  of  this  pair  of  joint  angles  will  not  change  the  result).  This 
result  is  applicable  to  the  entire  possible  range  of  the  two  Dundur's  parameters.  For  this  specific 
pair  of  joint  angles,  the  stress  singularity  is  limited  to  a  very  small  zone  near  a  =1.  These 
extreme  material  joint  combinations  are  quite  rare  in  engineering  applications  since  they 
represent  extremely  high  mismatch  in  Young’s  moduli. 

Planar  and  axisymmetric  bonded  specimens  (with  the  same  bonding  area)  were  used  for 
testing  specimens  in  tension.  The  planar  tension  specimens  are  254  mm  long  (individual  halves 
are  127  mm  long),  38  mm  wide  and  6  mm  thick.  These  specimens  were  used  to  make  same  and 
bi-material  joints.  The  length  and  thickness  of  all  the  specimens  were  chosen  such  that  the 
gripping  pressure  from  the  MTS  system  did  not  cause  any  specimen  damage.  The  two  types  of 
planar  tension  specimens  with  straight-edges  and  convex  edges  are  illustrated  in  Figure  3(a),  (b), 
and  the  two  types  of  axisymmetric  specimens  with  straight  edges  and  convex  edges  are  shown  in 
Figure  3(c),  (d).  The  free-edge  stress  singularity  was  completely  removed  in  axisymmetric 
convex  joint  specimens  (Wang  and  Xu,  2006). 
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Fig.3 


3.2  Comparisons  of  tensile  bonding  strengths  for  different  material  systems  and  adhesives 

In  each  of  the  tables  such  as  Table  1,  the  results  for  same-material  joints  are  reported  first 
followed  by  the  bi-material  data.  The  first  column  mentions  the  type  of  bonding  and  the  joint 
materials.  The  second  column  presents  the  data  as  a  mean  strength,  while  the  third  column 
reports  the  difference  between  the  weak  and  strong  bonding  of  the  same-material  joints.  In  all  the 
tables.  Loctite  adhesives  are  referred  to  by  their  corresponding  numbers.  For  example.  384  would 
denote  the  weak  bond  produced  by  utilizing  Loctite  384  adhesive.  W10  refers  to  the  strong 
adhesive  Weldon- 10.  This  convention  is  used  in  all  the  tables  presented  in  this  study. 

Tabic  1 

Material  systems  bonded  with  Weldon-10  show  a  higher  tensile  strength  than  the  same- 
material  systems  bonded  with  Loctite  384  as  seen  in  Table  1.  A  general  trend  observed  in  our 
tensile  strength  data  is  that  PMMA  material  systems  show  a  higher  value  of  tensile  strengths  than 
polycarbonate  systems.  All  of  the  specimens  were  observed  to  fail  in  sudden  and  brittle  tension 
along  the  interface.  To  ensure  repeatability,  at  least  five  specimens  were  tested  in  each  case.  In 
order  to  provide  a  complete  data  base  of  the  same-material  bonding  of  different  polymer 
systems,  some  results  from  our  previous  conference  paper  (Xu  et  ah,  2002)  are  presented  in 
Table  1  along  with  latest  data.  Experimental  data  from  our  previous  work  is  accumulated  here 
along  with  currently  measured  data  for  the  convenience  of  modeling  and  simulation  work. 

For  the  Homalite/adhesive/Homalite  system,  in  addition  to  the  strong  and  weak  adhesives 
(Weldon- 10  and  Loctite  384),  Polyester  was  also  used  in  strong  bonding.  The  tensile  strength  of 
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this  bond  was  the  highest  as  Polyester  shares  a  similar  chemical  structure  with  Homalite.  Loctite 
330  is  another  strong  adhesive  similar  to  Weldon- 10  and  its  tensile  bonding  strengths  are  quite 
high.  Indeed,  the  weak  adhesive  Loctite  384  only  provided  slightly  lower  bonding  strength  than 
Loctite  330  or  Weldon-10.  The  weakest  adhesive  used  is  Loctite  5083  since  its  tensile  strength  is 
only  20%  of  the  bonding  strength  obtained  by  using  the  strong  adhesive  Weldon- 10. 

Table  1  shows  a  slight  increase  in  the  tensile  bonding  strengths  of  bi-materials  for  those 
material  systems  using  Weldon- 10  adhesives.  Data  related  to  bi-material  bonding  are  from  our 
previous  measurements  (Wang  and  Xu,  2006).  Also,  bonding  strengths  of  the 
polycarbonate/aluminum  system  are  just  slightly  lower  than  those  of  PMMA/aluminum  system. 
However,  for  the  Homalite/steel  systems,  the  bonding  strength  using  Loctite  384  is  only  32%  of 
the  bonding  strength  of  PMMA/aluminum  system.  This  is  probably  due  to  the  possibility  that 
chemical  bonding  capability  of  Loctite  383  with  steel/Homalite  is  not  as  good  as  that  with 
aluminum/PMMA. 

4.  SHEAR  STRENGTHS  OF  SAME  AND  DISSIMILAR-MATERIAL  JOINTS 

4. 1  Specimens  and  test  approaches 

Shear  bonding  strength  is  a  key  material  property  for  any  adhesive  joints  of  dissimilar 
materials.  Previous  researchers  (Xu  and  Rosakis,  2002;  Kitey  and  Tippur.  2008)  have  shown  that 
shear  failure  is  a  dominant  interfacial  failure  mode  when  layered  materials  with  interfaces  were 
subjected  to  out-of-plane  dynamic  loading.  Thus,  for  dissimilar-material  joints  with  interfaces,  it 
becomes  very  necessary  to  measure  the  interfacial  shear  strength  for  quality  control  and  failure 
prevention.  In  our  experiments,  we  measured  the  interfacial  shear  strengths  over  a  range  of 
bonded  material  systems  using  two  types  of  shear  tests:  losipescu  and  short-beam  shear  tests 
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which  show  different  interfacial  shear  stress  distributions.  While  the  short-beam  shear  specimen 
shows  a  parabolic  variation  of  shear  stress  across  its  interface,  the  Iosipescu  specimen 
demonstrates  a  near  constant  shear  stress  variation  as  illustrated  in  Fig.4.  The  Iosipescu  set-up 
was  used  to  measure  the  shear  strength  of  both  types  of  specimens  w'ith  the  same  specimen 
width.  Iosipescu  and  short-beam  shear  specimens  were  76.2  mm  long  (individual  halves  were 
38.1  mm  long),  19.1  mm  wide  and  5.4  mm  thick.  In  addition,  Iosipescu  specimens  have  a  notch 
with  a  depth  of  3.8  mm  at  the  center  (Deng  et  al.,  2009).  Short-beam  shear  specimens  were 
exclusively  utilized  to  measure  the  shear  strength  of  the  bonded  same-materials.  Iosipescu  shear 
specimens  were  used  to  measure  the  shear  bonding  strength  for  both  same  and  bi-material  joints. 

Fig.4 

Coherent  Gradient  Sensing  experiments  were  carried  out  to  obtain  fringe  patterns  from 
bi-material  specimens  of  PMMA/W1 0/aluminum.  This  is  shown  in  Fig.  5  where  an  Iosipescu 
specimen  is  subjected  to  increasing  shear  load.  Aluminum  is  not  transparent  and  hence  the 
fringes  were  seen  only  on  the  PMMA  block.  The  set  of  four  pictures  shown  in  Fig.  5  illustrate 
the  fringe  development  in  the  specimen  from  a  small  initial  load  until  failure.  A  strong  stress 
concentration  is  seen  at  the  loading  points  and  near  the  concave  joint  between  aluminum  and 
PMMA.  Unlike  tensile  specimens,  we  were  unable  to  design  convex  joints  for  bi-material  shear 
specimens;  hence  the  free-edge  stress  singularity  still  existed  in  bi-material  shear  specimens. 

Fig.  5 

In  case  of  specimens  with  polycarbonate,  photoelasticity  experiments  were  carried  out  to 
obtain  fringe  patterns  which  are  contours  of  maximum  in-plane  shear  stress  given  by  equation 
(1).  These  patterns  as  seen  in  Fig.  6,  w’ere  obtained  for  short-beam  shear  specimens  at  loads  of 
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12.5  %  and  25  %  of  their  final  failure  load.  A  finite  element  model  of  the  specimen  was  built 
w  ith  AN  SYS  1 1.0  to  obtain  the  stress  distributions  in  the  bonded  specimen.  A  two  dimensional 
analysis  was  employed  for  bonded  polycarbonate  materials.  The  mesh  comprised  of  only  Plane 
42  elements  to  meet  the  requirements  of  the  plotting  software  Tecplot.  Displacements  were 
applied  to  the  specimen  edges  to  simulate  realistic  boundary  conditions.  An  iterative  procedure 
was  adopted  to  ensure  that  the  nodes  at  the  lower  boundary  were  free  of  tensile  load  as  it  is 
impossible  for  these  nodes  to  be  in  tension.  The  obtained  loading  pattern  was  anti-symmetric 
across  the  interface.  Further  details  of  the  loading  condition  and  the  finite  element  model  can  be 
found  in  Krishnan  and  Xu  (2009).  This  stress  value  was  then  converted  to  a  fringe  order  N  which 
in  turn  was  converted  into  a  grayscale  value.  Half  order  fringes  (0.5,  1 .5,  2.5,  etc.)  were  given  a 
value  of  255  and  full  order  fringes  (0,  1,  2.  etc.)  were  assigned  a  value  of  0  on  the  grayscale 
spectrum.  The  resulting  fringe  patterns,  plotted  with  Tecplot,  are  compared  with  the 
experimentally  observed  patterns  as  depicted  in  Fig.  6.  A  close  match  was  obtained  between  the 
simulated  and  experimental  patterns  thereby  validating  our  finite  element  model. 

Fig.  6 

4.2  Shear  strength  comparisons  and  effect  of  the  interface  shear  stress  on  the  bonding  strength. 

The  shear  strengths  from  the  different  material-systems  are  reported  in  Table  2  and  Table 
3.  Table  2  presents  shear  strength  data  obtained  exclusively  by  testing  losipescu  shear  specimens 
for  a  variety  of  indicated  material  systems.  In  each  case,  about  5-6  specimens  were  tested  to 
ensure  repeatability  and  the  mean  values  along  with  their  standard  deviations  are  reported  here. 
All  the  shear  specimens  were  found  to  fail  in  shear  along  the  bonding  line  in  a  sudden  and  brittle 
fashion  except  for  Homalite/WlO/Homalite  and  Homalite/polyestcr/Homalite  systems.  Bonded 
same-materials  consistently  show  a  greater  value  of  shear  strengths  in  comparison  with  bi- 
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material  systems  with  the  same  adhesive  bonding.  This  is  attributed  to  a  weaker  bonding  in  case 
of  bi-materials  due  to  the  mismatch  in  elastic  properties.  Similar  to  tensile  bonding  strengths, 
material  systems  with  PMMA  show  a  higher  value  of  mean  shear  strength  than  systems  with 
polycarbonate.  Again,  the  difference  between  Loctite  384  and  Weldon-10  bonding  is 
significantly  higher  for  PMMA  bonded  specimens  than  its  polycarbonate  counterparts  in  case  of 
the  same-material  bonding.  Materials  bonded  with  Weldon- 10  consistently  report  a  value  of 
higher  shear  strength  than  Loctite  384. 

Bonding  strengths  of  Homalite/adhesive/Homalite  systems  using  different  adhesives  are 
also  presented  in  Table  2.  Shear  bonding  strengths  of  the  same-material  systems  show  very 
different  values  for  various  adhesives  just  as  in  tensile  bonding  strengths  of  the  same 
Homalite/adhesive/Homalite  systems.  Homalite  shear  specimens  bonded  with  strong  adhesives 
Weldon- 10  and  polyester,  were  found  to  fail  in  tension  in  the  Homalite  part,  rather  than  in  shear 
across  the  bonded  interface  (Xu  et  al.,  2004b).  Thus,  its  actual  shear  bonding  strength  cannot  be 
measured;  but,  the  lower  limit  of  the  shear  bonding  strength  is  obtained  as  seen  in  Table  2. 
However,  for  PMMA  and  polycarbonate,  actual  shear  bonding  strengths  can  be  measured  using 
Iosipescu  shear  tests  as  their  tensile  strengths  are  much  higher  than  that  of  Homalite. 

Table  2. 

Table  3  contains  shear  strength  data  obtained  from  both  Iosipescu  and  short-beam  shear 
tests.  Here,  in  each  case  about  25-30  specimens  were  tested.  Also,  specimens  for  which  data  is 
reported  in  Table  3  were  bonded  only  using  weak  adhesive  Loctite  384.  These  data  indicate  only 
a  small  difference  between  measured  shear  bonding  strengths  of  Iosipescu  and  short-beam  shear 
specimens,  although  their  interfacial  shear  stress  distributions  are  very  different  as  shown  in 
Fig.4.  Therefore,  interfacial  shear  stress  distributions  have  the  least  effect  on  the  shear  strength 


12 


Krishnan  and  Xu  (2009) 


Bonding  strength  and  fracture  toughness  of  adhesive  joints 


measurement.  It  has  been  demonstrated  that  short-beam  shear  specimens  are  more  efficient  in 
machining  than  Iosipescu  shear  specimens  (Krishnan  and  Xu,  2009). 

Table  3. 


5.  FRACTURE  TOUGHNESS  FOR  SAME  AND  DISSIMILAR-MATERIAL  JOINTS 

5. 1  Principles  on  fracture  toughness  measurements 

Edge-notched  fracture  specimens  were  used  to  measure  the  mode  I  fracture  toughness. 
The  specimens  were  designed  and  tested  to  the  following  dimensions:  specimen  width  W=38 
mm,  specimen  thickness  B=5.4  mm  and  initial  crack  length  a=19  mm.  All  specimens  had  an 
initial  crack  with  aTT  =0.5  made  before  bonding  by  using  scotch  tape.  All  of  the  fracture 
specimens  were  tested  in  three-point  bending  (Suresh  et  al.,  1990).  The  length  of  the  initial  crack 
in  a  bi-material  specimen  in  such  a  case  would  have  to  be  chosen  carefully  to  enable  fracture  in  a 
single  mode  (mode  I).  The  mode  I  fracture  toughness  Kic  for  same-material  joints  was  calculated 
using  equation  (5)  (Anderson,  1995). 


K,c  = 


PQS 

BWvl 


f(x)  0<x  =  a/fV  <1 


fix) 


3  r-  [1 .99- jc(1  - jc)(2.15  -  3.93jr  +  2.7jr2)] 
2V*  (1  +  2jc)(1  -  jc)3/2 


(5) 


where  Pq  is  the  maximum  load  from  the  load-displacement  plot.  S’  is  the  support  span,  f(x) 
accounts  for  the  correction  due  to  the  specimen  geometry.  In  case  of  bi-materials,  the  calculation 
of  fracture  toughness  becomes  very  different.  First,  the  asymptotic  stress  field  of  an  interfacial 
crack  in  a  bi-material  specimen,  oy  can  be  expressed  as  (Rice,  1988) 


cr  = 


v2 


7T  V 


[Re{A'r"}<7  (6\s)+  Im  {  K  r  ,c }  cr  (0  ;  s  )] 


(6) 
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where  K=Ki+iK2  is  the  complex  stress  intensity  factor,  a,/  and  G,j11  are  the  stresses  in  mode  I  and 
mode  II,  and  e  is  a  function  of  Dundur’s  parameters  p  and  is  given  by 


s 


2  71 


In  { 


i  -  a 

i  +  p 


(7) 


The  elastic  properties  of  aluminum  include  a  Young’s  modulus  of  E=71  GPa,  a  shear  modulus 
p=26.7  GPa  and  a  Poisson's  ratio  of  v=0.33.  Corresponding  elastic  properties  for  polycarbonate 
used  in  this  calculation  are  E=2.4  GPa,  p=0.9  GPa,  v=0.34.  Hence,  the  two  Dundur’s  parameters 
for  the  material  combination  of  polycarbonate  and  aluminum  are  calculated  as  a=0.93  and 
P=0.31.  It  should  be  noted  here  that  PMMA,  although  chemically  different,  has  similar  elastic 
properties  as  polycarbonate.  Therefore,  the  Dundur's  parameters  for  PMMA/aluminum  are  also 
the  same  as  in  above.  A  schematic  of  our  bi-material  specimen  used  to  obtain  the  fracture 
toughness  value  is  illustrated  in  Fig.  7. 

Fig.7 

A  general  form  for  the  stress  intensity  factor  for  a  bi-material  specimen  is  given  as  (O'Dowd  et 
al..  1992) 


K  =  YT  4aa~ieeiv  (8) 

where  T=P(3BAV2),  Y  and  v| /  are  calibrating  factors  which  depend  on  a/W,  B/W and  the  Dundur’s 
parameters.  Then  the  stress  intensity  factor  in  mode  I  can  be  expressed  as 

K  ,  =  Re{KaiE}  (9) 

Using  equations  (8)  and  (9).  the  bi-material  fracture  toughness  Kic  is  calculated  as 
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KIC  =  —^~ry  x/tf  cos(^ )  0°) 

>c  BW2 

Here,  Y=2.4  and  vj/=7.8  degrees,  the  fitting  parameters  are  obtained  from  0‘Dowd  et  al.  (1992) 
for  P=ct/3  as  is  in  our  case. 

5.2  Fracture  toughness  results 

The  material  systems  and  their  measured  Kic  values  are  presented  in  Table  4.  During 
experiments,  the  load  required  to  completely  break  the  specimen  under  three-point  bending  was 
recorded  and  the  value  of  Kic  was  calculated  using  equation  (5)  for  same-material  joints  and 
using  equation  (10)  for  bi-material  systems.  It  was  observed  that  Weldon-10  bonding  shows  a 
slightly  higher  value  of  Kic  than  Loctite  384  bonding  for  most  material  systems.  Only 
polycarbonate/384/aluminum  shows  a  decrease  of  10%  in  fracture  toughness  from  the 
corresponding  weak  bond.  However,  it  was  also  observed  that  the  standard  deviation  in  case  of 
Weldon- 10  bonding  was  higher  than  Loctite-384  bonding.  Generally  bi-material  fracture 
specimens  show  a  small  difference  of  Kic  values  between  the  two  types  of  strong  and  weak 
adhesives.  Similar  to  its  tensile  and  shear  bonding  strengths.  PMMA  shows  a  better  bonding  with 
the  two  ty  pes  of  adhesives  and  hence  a  greater  value  of  fracture  toughness  is  obtained  for  bonded 
PMMA  specimens.  Bi-material  specimens  consistently  show  a  lower  value  of  mode  1  fracture 
toughness  than  same-material  specimens. 

Surprisingly,  Homalite/polyester/Homalite  material  systems  show  lower  fracture 
toughness  than  other  strong  adhesive  systems  used  in  conjunction  with  Homalite:  Weldon-10  and 
Loctite  330.  It  should  be  noted  that  polyester  provides  the  highest  tensile  and  shear  bonding 
strengths  as  can  be  seen  in  Tables  1  and  2.  Therefore,  strength  and  fracture  toughness  are  very 
different  parameters  and  should  be  measured  for  every  new  material  system.  The  fracture 
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toughnesses  of  other  adhesive  systems  show  a  trend  similar  to  the  bonding  strengths  of  the  same- 
material  joints,  i.e.,  bonding  strengths  increase  with  fracture  toughnesses  from  weak  bonding  to 
strong  bonding. 

Table  4 

6.  CONCLUSIONS 

Five  different  types  of  strong  and  weak  adhesives  (Weldon- 10,  Loctite  384,  Loctite  330, 
Loctite  5083  and  polyester)  were  used  in  conjunction  with  five  types  of  materials  (aluminum, 
steel,  PMMA,  polycarbonate  and  Homalite)  to  produce  a  variety  of  bonded  material  systems. 
These  material  systems  (same  and  bi-material  joint)  are  tested  in  shear,  tension  and  fracture. 
Results  indicate  that  materials  bonded  with  Weldon-10  and  polyester  consistently  show  higher 
tensile  and  shear  bonding  strengths  than  the  same-material  systems  bonded  with  other  adhesives. 
However,  bi-material  fracture  toughness  (mode  I)  tends  to  show  lesser  difference  between  the 
weak  and  strong  adhesive  bonding.  In  general,  bi-material  systems  in  shear  and  fracture  report 
lower  properties  than  the  same-material  systems  due  to  higher  property  mismatch  involved  in  bi¬ 
material  systems.  Specimens  bonded  with  PMMA  consistently  show  higher  values  of  shear  and 
tensile  bonding  strengths,  and  mode  1  fracture  toughness  than  polycarbonate  systems.  A  wide 
variety  of  carefully  measured  experimental  data  will  be  highly  beneficial  for  computational 
simulations  on  failure  mechanics  of  adhesive  joints. 
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Figure  1.  Experimental  setup  of  a  mechanical-optical  system  to  record  in-situ 
photoelasticity  during  loading  process 
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Figure  2.  (a)  Coherent  Gradient  Sensing  (CGS)  photograph  showing  a  strong  stress 
concentration  (associated  with  fringe  pattern  concentrations)  at  the  free  edges  of  a  bonded 
metal  and  polymer  subjected  to  tensile  load  (b)  Angular  definition  of  a  bi-material 
wedge. 
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Figure  3.  Schematic  diagrams  of  metal-polymer  joint  specimens  with  (a)  straight 
edges  (baseline);  (b)  convex  edges  with  least  stress  singularities;  (c)  axisymmetric 
straight  joints  (baseline);  (d)  axisymmetric  convex  joints  with  least  stress  singularities. 
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Figure  4.  Schematic  models  indicating  the  structure  of  (a)  Short-beam  shear  and  (b) 
Iosipescu  specimen. 
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Figure  5.  CGS  fringe  pattern  development  for  bonded  PMMA/W  10/aluminum  losipescu 
shear  specimen.  The  First  snapshot  shows  fringe  just  evolving  as  the  loading  is  increased 
and  the  last  snapshot  shows  the  final  fracture  in  the  bonded  bi-material  specimen  with  a 
clear  shear  displacement  at  the  interface. 
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Figure  6.  Comparison  between  experimental  and  finite  element  generated  fringe  patterns 
for  a  typical  polycarbonate  short-beam  shear  specimen. 
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Figure  7.  Bi-material  specimen  for  mode-1  fracture  toughness  measurement  (a/W=0.5) 
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Table  1.  Tensile  bonding  strengths  for  same  and  bi-material  joints 


Material/adhesive/material 

Tensile  Strength  (MPa) 

Notes 

Same-Material  Bond 

Polycarbonate  /384/  Polycarbonate 

6.06  ±  1.64 

+  113%  increase 
strong/weak 

Polycarbonate  AV 1 0/  Polycarbonate 

1 2.93  ±  1.36 

PMMA/384/PMMA 

12.66  ±  1.63 

+64  %  increase 
strong/weak 

PMMAAV 1 0/PMMA 

20.87  ±6.21 

Homalite/Polyester/Homalite 

28 

Xu  et  al.,  2002 

HomaliteAVl  0/Homalite 

7.74 

Homalite/330/Homalite 

7.0 

Homalite/384/Homalite 

6.75 

Homalite/5083/Homalite 

1.53 

Xu  et  al.,  2002 

Bi-Material  Bond 

Polycarbonate  /384/aluminum 

9.57  ±0.72 

Polycarbonate  /W  10/aluminum 

11.35 

Wang  and  Xu, 
2006 

PMMA/3  84/aluminum 

10  01  ±0.93 

Wang  and  Xu, 

PMMAAV  10/  aluminum 

12.85 

2006 

Homalite/330/steel 

5.38 

Xu  et  al.,  2002 

Homalite/384/steel 

3.25 
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Table  2.  Shear  bonding  strengths  for  same  and  bi-material  joints 


Material/adhesive/material 

Shear  Strength  (MPa) 

Notes 

Same-Material  Bond 

Polycarbonate/384/  Polycarbonate 

10.99  ±  1.45 

+41  % 

Po ly carbonate AV 10/  Poly  carbonate 

15.52  ±2.74 

PMM  A/3  84/PMM  A 

1 1.58  ±  2.15 

+  119% 

PMMA/W 1 0/PMMA 

25.35  ±  2.25 

Homalite/WlO/Homalite 

>21.65 

Xu  et  al.,  2002 

HomaliteTolyester/Homalite 

>23.26 

Homalite/330/Homalite 

12.58 

Homal  ite/384/Homalite 

7.47 

Homalite/5083/Homalite 

0.81 

Xu  et  al.,  2002 

Bi-Material  Bond 

Polycarbonate/3  84/al  um  i  num 

6.18  ±  2.16 

65  % 

PolycarbonateAV  10/aluminum 

10.1 8  ±  4.35 

PMMA/384/aluminum 

8.63  ±  0  85 

18% 

PMMA/WD  10/aluminum 

10  16  ±3.92 
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Table  3.  Shear  bonding  strengths  of  Iosipescu  and  short-beam  shear  specimens 


Material/adhesive/material 

Iosipescu  Shear  (MPa) 

Short-Beam  Shear(MPa) 

Difference 

Aluminum/3  84/Aluminum 

10.75  ±2.39 

10.16  ±  2.41 

5.5% 

Polycarbonate/3  84/Poly  carbonate 

10.99  ±  1.45 

8.51  ±  1.13 

22.5  % 

PMMA7384'PMMA 

1 1.58  ±  2.15 

10.19  ±  0.57 

12% 
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Table  4.  Mode-1  fracture  toughness  for  same  and  bi-material  joints 


Material  adhesive/material 

Mean  Klc(MPam17) 

Difference  (%) 

Same-Material  Bond 

Polycarbonate/3  84/Polycarbonate 

0.64  ±0.02 

±36  % 

Polycarbonate/W  1 0/Polycarbonate 

0.86  ±0.16 

PMMA/384/PMMA 

0.71  ±0.06 

±147% 

PMM  A/W 1 0/PMMA 

1.74  ±  0.15 

Homalite/Polyester/Homalite 

0.56 

Xu  et  al.,  2002 

Homalite/WlO/Homalite 

0.83 

Homalite/330/Homalite 

0.75 

Homalite/384/Homalite 

0.38 

Homalite/5083/Homalite 

0.10 

Xu  et  al.,  2002 

Bi-Material  Bond 

Polycarbonate/384/Aluminum 

0.1  ±0.01 

-10% 

Po  ly  carbon  ate/W  1 0/Aluminum 

0.09  ±0.02 

PMMA/3  84/Aluminum 

0  12  ±0.01 

±17% 

PMMA/W 1 0/  A 1  u  m  i  n  u  m 

0.14  ±0.02 
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